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1 Introduction 

Chemistry is primarily concerned not with the properties of single molecules but 
with periodic trends, homologous series and the like. It is, therefore, important 
that any method which we apply to the problem of molecular electronic structure 
depends linearly on the number of electrons in the system being studied. Mean- 
ingful comparisons of atoms and molecules of different sizes are then possible. 
This property has been termed size-consistency 1 - 2 . Independent electron models, 
such as the widely used Hartree-Fock approximation, provide a size-consistent 
theory of atomic and molecular structure. 

Independent-electron models account for the major proportion, typically 
99.5%, of the non-relativistic electronic energy of an atom or molecule. The 
Hartree-Fock model describes not only the Fermi interactions of the electrons but 
also their averaged electrostatic interactions. It is unfortunate that the remaining 
energy is of the same order of magnitude as most energies of chemical interest. 
This remaining energy, the correlation energy, arises from the 'instantaneous 
correlations' of the individual electronic motions. Chemistry is primarily 
concerned with small energy differences, such as those between different nuclear 
geometries or different electronic states, and these differences may be seriously 
affected by the correlation energy. 

In the past twenty years, there has been increasing interest in the calculation 
of correlation energies and other properties of atomic and molecular systems by 
means of diagrammatic many-body perturbation theory techniques 3-9 due to 
Brueckner 10 and Goldstone. 11 Diagrammatic many-body perturbation theory 
provides a simple pictorial representation of electron correlation effects in atoms 
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Mechanics', Cambridge University Press, 1967. 

4 A. L. Fetter and J. D. Walecka, 'Quantum Theory of Many-particle Systems', McGraw- 
Hill, New York, 1971. 

5 H. P. Kelly, Adv. Chem. Phys., 1969, 14, 129. 

6 J. Paldus and J. Cizek, Adv. Quantum Chem., 1975, 9, 105. 

7 I. Hubac and P. Carsky, Top. Curr. Chem., 1978, 75, 97. 

8 S. Wilson, in 'Proceedings of Daresbury Study Weekend', Dec. 1977, ed. V. R. Saunders, 
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and molecules and also forms the basis of a tractable, non-iterative scheme for 
accurate calculations. 5 - 12-14 Perturbation theory provides perhaps the most 
systematic technique for the evaluation of corrections to independent electron 
models. 

The many-body perturbation theory is so called because it can be applied to 
arbitrarily large systems. In fact, the theory was originally devised to treat infinite 
fermion systems. It leads to expressions for correlation corrections to indepen- 
dent electron models which have a linear dependence on the number of electrons 
being considered. If the theory is applied to a system A, giving an energy E(A), 
and to a system B, giving an energy E(B), and then to the combined system AB, 
where A and B are an infinite distance apart, then the energy of the super- 
system, £XAB), is given by , 

E(AB) = E(A) + E(B) (1) 

The energy of any system may be written as a sum of the energies of its component 
parts no matter how these components are defined. This property is not shared by 
some of the other methods currently employed in the study of electron correla- 
tion, for example the widely used method of configuration mixing limited to 
single- and double-excitations, which, when a single determinantal reference 
function is used, leads to an expression for the correlation energy depending on 
the square root of the number of electrons under consideration. 15 - 16 Limited 
configuration mixing is not a size-consistent technique. 

The diagrammatic many-body perturbation theory may be derived from the 
Rayleigh-Schrodinger perturbation expansion. Brueckner 10 showed that certain 
terms arise in the Rayleigh-Schrodinger expansion which have a non-linear 
dependence on the number of electrons being studied. He showed that these 
unphysical terms cancel in each of the first few orders of the Rayleigh-Schrodin- 
ger perturbation series. Goldstone 11 generalized this result to all orders using 
the diagrammatic techniques of time-dependent perturbation theory. This leads 
to the linked diagram perturbation expansion. 3-7 All terms corresponding to 
unlinked diagrams depend non-linearly on the number of electrons and thus 
mutually cancel in each order. This cancellation of unlinked diagrams not only 
eliminates unphysical terms but also leads to important computational simplica- 
cations. 12-14 

The pioneering work on the application of the many-body perturbation 
theory to atomic and molecular systems was performed by Kelly. 5 - 17-21 He 
applied the method to atoms using numerical solutions of the Hartree-Fock 
equations. Many other calculations on atomic systems were subsequently 
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80 H. P. Kelly, Adv. Theor. Phys., 1968, 2, 75. 
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reported (e.g. refs. 22 — 26). The first molecular calculations using many-body 
perturbation theory used single-centre expansions and were limited to simple 
hydrides where it is possible to treat the hydrogen atoms as additional perturba- 
tions. 27-31 More recently, the theory has been applied to arbitrary molecules by 
employing the algebraic approximation 32-42 which is fundamental to most 
molecular calculations. In this approximation, single-particle state functions are 
parameterized in terms of a finite basis set. This is equivalent to replacing the 
true hamiltonian by a model hamiltonian whose domain is restricted to some 
subspace of the Hilbert space associated with the true hamiltonian. 

In this article, the results of atomic calculations will only be considered when 
they are relevant to the molecular situation. This is the case in a number of areas 
where the application to atoms is well established but remains to be extended 
to molecules. This article is concerned with the application of the many-body 
perturbation theory to arbitrary molecular systems. Recent work 43-44 has shown 
that this technique can be at least as if not more accurate than other techniques 
currently employed in the study of molecular electronic structure. The method 
is probably computationally more efficient than other schemes and certainly 
has a number of theoretical properties which make its use attractive. For example, 
in discussing the widely used method of configuration mixing, Shavitt states : 46 
'The fact that in a configuration interaction expansion unlinked cluster contribu- 
tions can only be accounted for by including quadruple- (and higher-order) 
excitations is one of the principal drawbacks of the method. In contrast, such 
contributions are automatically accounted for without explicitly computing 
higher-order terms, in some cluster-based methods and in many-body perturba- 
tion theory. In this sense the CI expansion is much less compact and less efficient 
than these approaches and becomes progressively less efficient as the number of 
electrons increases.' 

22 E. S. Chang, R. T. Pu, and T. P. Das, Phys. Rev., 1968, 174, I. 

23 N. C. Dutta, C. Matsubara, R. T. Pu, and T. P. Das, Phys. Rev., 1969, 177, 33. 

24 R. T. Pu and E. S. Chang, Phys. Rev., 1966, 151, 31. 

25 T. Lee, N. C. Dutta, and T. P. Das, Phys. Rev. A., 1970, 1, 995. 

26 T. Lee, N. C. Dutta, and T. P. Das, Phys. Rev. A., 1971, 4, 1410. 

27 H. P. Kelly, Phys. Rev. Lett., 1969, 23, 455. 

28 J. H. Miller and H. P. Kelly, Phys. Rev. Lett., 1971, 26, 679. 

29 T. Lee, N. C. Dutta, and T. P. Das, Phys. Rev. Lett., 1970, 25, 204. 

30 T. Lee and T. P. Das, Phys. Rev. A., 1972, 6, 968. 

31 C. M. Dutta, N. C. Dutta, and T. P. Das, Phys. Rev. Lett., 1970, 25, 1695. 

32 S. Wilson and D. M. Silver, Phys. Rev. A., 1976, 14, 1949. 

33 J. M. Schulman and D. N. Kaufman, /. Chem. Phys., 1970, 53, 477. 

34 J. M. Schulman and D. N. Kaufman, /. Chem. Phys., 1972, 57, 2328. 
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37 M. A. Robb, Chem. Phys. Lett., 1973, 20, 274; and in 'Computational Techniques in 
Quantum Chemistry and Molecular Physics', ed. G. H. F. Diercksen, B. T. Sutcliffe, and 
A. Veillard, D. Reidel, 1974, p. 435. 
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8 » R. J. Bartlett and D. M. Silver, /. Chem. Phys., 1975, 62, 3258; erratum, 1976, 64, 4578. 

40 D. F. Freeman and M. Karplus, /. Chem. Phys., 1976, 64, 2641. 

41 M. Urban, V. Kello, and I. Huba5, Chem. Phys. Lett., 1977, 51, 170. 

42 S. Prime and M. A. Robb, Theor. Chim. Acta, 1976, 42, 181. 

43 S. Wilson and D. M. Silver, /. Chem. Phys., 1977, 66, 5400. 

44 S. Wilson and D. M. Silver, /. Chem. Phys., 1977, 67, 1649. 
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Theory', ed. H. F. Schaefer III, Plenum Press, New York, 1977. 
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This article is divided into seven parts. The many-body perturbation theory is 
discussed in the next section. The algebraic approximation is discussed in some 
detail in section 3 since this approximation is fundamental to most molecular 
applications. In the fourth section, the truncation of the many-body perturbation 
series is discussed, and, since other approaches to the many-electron correlation 
problem may be regarded as different ways of truncating the many-body pertur- 
bation expansion, we briefly discuss the relation to other approaches. Computa- 
tional aspects of many-body perturbative calculations are considered in section 5. 
In section 6, some typical applications to molecules are given. In the final section, 
some other aspects of the many-body perturbation theory of molecules are 
briefly discussed and possible directions for future investigations are outlined. 

2 The Many-body Perturbation Theory 

General Remarks. — In this section a brief introduction to the many-body per- 
turbation theory is given. In the second part the partitioning technique due to 
Lowdin 4 * and Feshbach 47 is used to give a straightforward and general introduc- 
tion to perturbation expansions. The perturbation series of Lennard-Jones, 48 
Brillouin, 49 and Wigner 50 is then described. This series is not suitable for applica- 
tion to many-particle systems and we, therefore, indicate how the many-body 
perturbation theory can be derived from the Rayleigh-Schrodinger perturbation 
theory. Diagrammatic rules and conventions are then introduced enabling the 
diagrammatic formulation of the many-body perturbation theory to be given. 
Some generalizations of the theory are briefly considered in the final part of this 
section. 

The Partitioning Technique. — Let P denote the projector onto some zero-order 
model wave function | <£„> and Q its complement. The electronic Schrodinger 
equation 

3$ \wy = § |f> (2) 
may then be written as a two by two block matrix equation 
IP3&P Pj?Q\ (P r F\ _ (PV\ 

where|^o>=P|y / >-Q|y / > can now be eliminated to produce the effective 
Schrodinger equation 

[PJ?P + PJ&Q (<?o - QJPQY 1 QJ&P] |$o> = <?o| #o> (4) 

or 

3? e tt\&o> = <?o|#o> (5) 

where 

^eff = [PsfP + PrfQ (<?o- QrfQY 1 QJ&P] (6) 

46 P. O. Lowdin, J. Math. Phys., 1962, 3, 969 and references therein. 

47 H. Feshbach, Ann. Phys. (N. Y.), 1962, 19, 287. 

48 J. E. Lennard-Jones, Proc. R. Soc. London, Ser. A., 1930, 129, 598. 

49 L. Brillouin, /. Physique, 1932, 7, 373. 

50 E. P. Wigner, Math. u. naturw. Anz. ungar. Akad. Wiss., 1935, 53, 475. 
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This effective hamiltonian has eigenfunctions in the model space but has the 
exact energy as an eigenvalue. 

Various forms of perturbation theory result from different expansions of the 
inverse in the effective hamiltonian using the identity 

, n=oo 

{A-B)~ x = S A^iBA^Y (7) 
*=0 

If <# denotes some zero-order hamiltonian and E its ground state eigenvalue 
then the perturbation series of Lennard-Jones, 48 Brillouin, 49 and Wigner 50 is 
obtained by putting 

A = ^o-^o (8) 

and 

B = jT-j^o (9) 
The Rayleigh-Schrodinger perturbation expansion is obtained by putting 

A = Eo-^o (10) 

and 

B = 3#-3#o-<? + Eo (11) 

Lennard-Jones Brillouin Wigner Perturbation Theory. — Let us write the total 
hamiltonian operator as a sum of a zero-order operator and a perturbation 

= jp + 3& x (12) 

with 

■**o|<Pi> = Ei\^i> (13) 

and 

jt\v t > = <?,|^> = (Ei + AEkmy (14) 
We introduce the projection operators 

p = |<P ><4>o| (15) 

and 

Go = I-Po (16) 
and employ the intermediate normalization convention 

<^o|^o> = 1 = <#o|#o> (17) 
Now we can define the wave operator, with the following properties 

|^o> = 3|<Po> (18) 

and 

Po£ = Po (19) 
8Po = & (20) 
GQo = (21) 
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and thus obtain an expression for the level shift, AE , for the ground state 

/= 00 

i=l 

= <«Po|^i|^o> (22) 

= <&o\3?lG\<Po> 
= <0 O |^|#O> 

where the reaction operator, i^, is given by 

= 3#x& (23) 

In Lennard-Jones Brillouin Wigner perturbation theory the wave operator is 
written as 

and the level shift has the form 

A£„ = M [, + (—^-yrf,)"] |*„> (25) 

These expressions for the wave operator and the reaction operator are formally 
equivalent to the integral equations 



and 



a - p ° + £o + aV^ (26) 



+ (27) 



from which the corresponding perturbation expansions can be obtained by 
iteration. 

Explicitly, the first few terms in the Lennard-Jones Brillouin Wigner perturba- 
tion series take the form 

Eo = <$o|^o|$o> (28a) 
E< 1 > = <<*>o|^i|#o> (28b) 



£<2) = <p o \jtr l aue l \0o> (28c) 

£< 3 > = <<Po |^i(^^i) 2 1 0o> (28d) 

= <«Po |^i(^i) 3 1 <P > (28e) 

E<j» = <<P l^i(^^i) 4 1 #o> (28f) 



where & is the resolvent 

Go 



(29) 



The Lennard-Jones Brillouin Wigner perturbation expansion is a simple geo- 
metric series. However, it contains the unknown exact energy within the denomi- 
nators. This expansion is, therefore, not a simple power series in the perturbation. 
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The perturbation theory of Lennard- Jones, Brillouin, and Wigner is not size 
consistent. 

Rayleigh-Schrodinger Perturbation Theory. — In Rayleigh-Schrodinger perturba- 
tion theory the unknown energy in the denominators of the Lennard-Jones 
Brillouin Wigner expansion is avoided. This enables a size-consistent theory to 
be derived. 

The wave operator, i5, may be written in an alternative form by replacing 
byJ? + Q &E Q and^i by A - Q Q &E Q giving 

& = Po + i- A£b) Q (30) 

Rearranging this expression in terms of powers of the perturbation, we 
obtain 

The y-th order energy is given by 

£<» = <0o\jPiAI*-»\0o> (32) 

Explicitly, the first few terms in the Rayleigh-Schrodinger perturbation expansion 
may be written in the form: 

E = <C>o|^o|<Po> (33a) 

£<i> = <0>o|^i|#o> (33b) 

£<2> = <3> 1 j^i^o^i | (P > (33c) 

£(3) = <<P 1 ^(^o^i) 8 1 *o> - ^K^Po | J^i^gJ^i | *o> (33d) 
£(«) = <<P [^i(^o^i) 3 |^o>-£: a> {<*o|^i^o^i^^i|^> 

+ <*o|^i^g^i^o^i|*o>}-£ ( 2) <«Po|^i^o^i|<Po> (33e) 



where $ is the reduced resolvent 

Go 



(34) 



Clearly, the terms other than the first in the expressions for E ( l\ £<*>, E ( l\ etc. 
depend on the number of electrons in a non-linear fashion. These terms exactly 
cancel components of the first terms in each of the expressions which also have a 
non-linear dependence on the number of electrons. 

The Many-body Perturbation Theory. — The Rayleigh-Schrodinger form of per- 
turbation theory provides an expansion for expectation values which have a 
linear dependence on the number of electrons in the system, N. In each order, 
other than zero-, first-, and second-order, terms arise which have a non-linear 
dependence on N. Brueckner 10 showed that for the first few orders the terms 
having a non-linear dependence on N mutually cancel in each order. Goldstone 11 
showed, using time-dependent perturbation theory, that this result can be 
generalized to all orders. The terms having a non-linear dependence on N may 
be associated with unlinked diagrams while those having the desired linear 
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dependence on N are associated with linked diagrams. This is the well known 
linked diagram theorem of many-body perturbation theory. 

It should perhaps be stated at this point that the use of diagrams in the many- 
body perturbation theory is not obligatory. The whole of the theoretical apparatus 
can be set up in entirely algebraic terms. However, the diagrams are both more 
physical and easier to handle than the algebraic expressions and it is well worth 
the effort required to familiarize oneself with the diagrammatic rules and con- 
ventions. 

The linked diagram expansion has, indeed, been derived by many authors and 
we shall, therefore, content ourselves with a brief outline of the Goldstone 
derivation here referring the interested reader elsewhere for full details. 3-5 

Before outlining the Goldstone treatment, we shall briefly mention some other 
derivations of the linked diagram theorem. Of particular interest is the derivation 
given by Brandow 51 which is based on the expansion of the energy-dependent 
denominators in the Lennard-Jones Brillouin Wigner perturbation theory. 
Paldus and Cizek 6 have given a time-independent derivation using a generaliza- 
tion of Wick's theorem 62 for time-independent problems. This approach has also 
been followed by Hubac and Carsky. 7 

The many-body perturbation theory is developed in terms of some set of single 
particle states, <f> Pi which are eigenfunctions of some single-particle operator, /, 

/fyp = ep^p (35) 

with eigenvalues e p . In the second-quantized formalism the zero-order hamil- 
tonian has the form 

= j dn + + (n)/(n) +(n) (36) 

and the perturbation operator may be written as 

#\ = i j dn j dr 2 ++(n) * + (r 2 )g(r u r 2 ) 4<r2) «K' - j dn ++(n) V(n) +(r l) (37) 

where v + ('"i) and v'(' , i) are the usual creation and annihilation field operators, 
g(ri> »*2) is the two-electron potential, and V(r x ) is the effective potential which is 
added to the bare-nucleus hamiltonian to give the one-electron operator /(i^). 
There is, of course, considerable freedom in the choice of the effective potential. 

Use of the interaction representation in time-dependent perturbation theory 
and an adiabatic switching, (|a|f), of the perturbation yields the evolution 
operator 

n = oo 

0*{t, -oo) = /+ s 0fT(/,-co) (38) 

w=l 

where / is the identity operator and 0%(t,— oo) is proportional to the wth power 
of the perturbation: 

OS (t, - oo) = ( - i)» f ' d/i f tl d/ a . . . P"" 1 d/„ J^l (tdsPi (/a) . . . J&X (tn) 

J — 00 J — 00 J — 00 

(39) 

51 B. H. Brandow, Rev. Mod. Phys., 1967, 39, 771. 
G. C. Wick, Phys. Rev., 1950, 80, 268. 
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= exp (iJ^ /)^i exp (\<x\t) exp (-iJfV) (40) 



The function 



|y ° > = lis; <^i^ ft -col *o> L (41) 



which obeys the intermediate normalization condition <?P | < F > = 1 is a well- 
defined eigenfunction of the total hamiltonian. 53 The total energy is given by 

E=Eo + AE (42) 

where the level shift is 

L a _ <M$ a (t, -00)|<P > J t-0 

The denominator in this equation may be written as 

exp [<<P | 0" (0, - oo) | <Z\>>]l (44) 

where the subscript L indicates that only terms described by linked diagrams are 
to be included. From equations (41) and (44), we obtain after taking the limit, the 
Goldstone expression for the wave function 

|fo> = (0(O,-oo)|<Po»l (45) 

and the level shift 

A£ = <<Po 1^1(0)^(0, - oo) | #o>l (46) 

Analysis of the products of field operators in these equations leads to a 
representation of the wave function and of the level shift in terms of diagrams 
of the type first introduced by Feynman. These diagrams provide a simple 
pictorial description of electron correlation effects in terms of the particle-hole 
formalism. 

Diagrammatic Conventions. — A diagram is a device through which a correspond- 
ing algebraic expression can be obtained. A number of different diagrammatic 
conventions are in use. We shall follow the work of Brandow 51 below and then 
discuss briefly its relation with other commonly used conventions. 

The basic elements of the diagrams are shown in Figure 1. Figure 1 (a) shows 
the diagrammatic representation of a one-electron operator matrix element. 
Figure 1 (b) shows the representation of a two-electron matrix which in the 
Brandow scheme includes permutation of the two electrons involved. Upward 
(downward) directed lines represent particles (holes) created above (below) the 
Fermi level when an electron is excited. 

The rules for performing the translation from a given diagram of the Brandow 
form to the corresponding algebraic expression are as follows: 

(i) Label each downward directed line with a unique 'hole' index: ij, k, . . . 
and label each upward directed line with a unique 'particle' index a, b, c, 

(ii) There is a summation over each unique hole and particle index covering 
all permissible values of these indices. 



■» M. Gell-mann and F. Low, Phys. Rev., 1951, 84, 350. 
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Figure 1 Diagram elements: (a) one-electron operator, (b) two-electron operator, (c) par- 
ticle line, (d) hole line 

(iii) The numerator of the summand consists of a product of antisymmetrized 
two-electron integrals 

Sp Q r, = j dn j dr 2 ^(n) ^(r 2 ) n* {Mn) U^) - Un) Mn)} (47) 

where ^ is a one-electron function and r 12 is the inter-electronic separation. For 
each interaction, i.e. horizontal dashed line, there is an integral of this type. The 
orbital indices can be read from the labelled diagram. Hence the indices p, q, r, s, 
should correspond to hole or particle lines entering or leaving the interaction in 
the following order: 'left-out, right-out, left-in, right-in', respectively. 

(iv) The denominator of the summand consists of a product of factors. There 
are «— 1 factors corresponding to an nth order diagram such that there is a 
denominator factor arising between adjacent interactions in the diagram. A 
denominator factor consists of the terms 

SfrSe, (48) 

i a 

where the first summation is over all hole lines that extend between the adjacent 
interactions and the second summation is over all particle lines that extend 
between the interactions. The e are the one-electron orbital energies. 

(v) There is a multiplicative factor of 1/2 for each pair of 'equivalent' lines. 
An equivalent pair of lines is defined to be two lines beginning at one interaction 
and ending at another and both going in the same direction. 

(vi) There is a multiplicative factor of (— l) p , where p is the sum of h and /. 
h is the number of unique hole lines in the diagram and / is the number of fermion 
loops. A fermion loop is determined by following the hole and particle lines in 
the direction of the arrows to form a continuous closed loop. 
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As an example of the application of these rules consider the diagram represent- 
ing the third-order 'hole-hole' energy, E 3 (h-h), shown in Figure 2(a). The lines 
in this diagram may be labelled to give Figure 2(b). The lowest interaction line 
gives rise to the integral /yab, the middle line implies the integral /mj, and the 
upper line the integral fabki. The lower half of the diagram yields the denomi- 
nator factor ei+ej—ea—eb, while the upper half corresponds to a denominator 





Figure 2 Diagrammatic representation of the third-order ' hole-particle' energy 

factor of £*+£!— e a — e&. The lines labelled i and j, k and /, and a and b are 
equivalent pairs giving a factor of (l/2) 3 = 1/8. There are two fermion loops and 
four unique hole lines giving a factor of (- 1) 6 = 1. E 3 (h-h) thus corresponds to 
the algebraic expression 

l V V <f i}ab# klij^abkl fA( n 

where 

®pq...T8... = e p + e q + . . . — e a — e&— . . . (50) 

Two other diagrammatic conventions are commonly used in perturbation 
theory. The Goldstone diagrams 11 are similar to those of Brandow 61 except that 
the interaction lines do not include permutation of the two electrons involved. 
Thus there is a set of Goldstone diagrams, which are related by electron exchange, 
corresponding to each Brandow diagram. The diagrams of Hugenholtz 64 are in 
one-to-one correspondence with those of Brandow. The Hugenholtz diagrams 
can be obtained from the Brandow diagrams by replacing the interaction lines in 
the latter by a single dot. This correspondence is illustrated in Figure 3. 






Figure 3 Correspondence between Brandow and Hugenholtz diagrams 
" N. H. Hugenholtz, Physica, 1957, 23, 481. 
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Diagrammatic Perturbation Theory. — The zero-order and first-order terms in the 
diagrammatic perturbation theory for a closed-shell molecule, which is described 
by a single determinantal wave function in zero-order, sum to the expectation 
value of the total hamiltonian for that wave function. Thus for a Hartree-Fock 
reference function the sum of the zero-order and the first-order energies is the 
self-consistent-field energy. The correlation energy is given by the sum of the 
second-order and higher-order terms. 

The second-order and third-order diagrams, using the Brandow antisym- 
metrized vertices convention, are displayed in Figure 4. When the Hartree-Fock 
reference function is employed, there is one second-order diagram and three 
third-order diagrams, which are distinguished by the type of central interaction 
line : the particle-particle (p-p) diagram, the hole-particle (h-p) diagram, and the 
hole-hole (h-h) diagram. The second-order and third-order (p-p) diagrams have 
only two hole lines; the correlation effects which they describe are, therefore, 
exclusively two-body. The third-order (h-p) diagram contains three hole lines 
and can therefore describe three-body interactions; however, if two of the hole 
lines are associated with the same single-particle state, then it describes two-body 
correlations. Finally, the third-order (h-h) diagram has four hole lines and can 
describe two-body, three-body, and four-body interactions. 




E 2 E 3 (pp) E 3 (hp) E 3 (hh) 

Figure 4 Second- and third-order diagrams which arise when the Hartree-Fock model is 
used to obtain a reference function 

Explicit expressions corresponding to the four diagrams shown in Figure 4 can 
be written down by following the rules given in the previous section. The second- 
order energy expression has the form 

E 2 = is S t mbf T (51) 

ij ab ^tjab — dijab 

where the dyab will be discussed below but in the immediate discussion should be 
taken to be zero. The third-order (p-p) energy may be written 

Ez(p-p) = I S S -iw ^ , ^ fa Zt fcai> l i ^ (52) 

ij abed (y?i]al>-di}abX@Vcd — di]cd) v 

the third-order (h-p) energy written as 

E 3 (h-p) = S S fS) fiia } f ^% f - bcile , v (53) 

Ok abc \&Uab — dijabX&)kbc — ajkbe) 

and the third-order (h-h) energy written as 

£ 3 (h-h) = | S S75^4l^V- f (54) 

tiki ab K&im-di)ab)(j20klab — dklab) K ' 
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Note that all of the above expressions are written in terms of single electron 
functions and no reference is made to many-electron functions. This is a funda- 
mental characteristic of the many-body perturbation theoretic approach to the 
correlation problem. 

In the above expressions the JV-electron Hartree-Fock model hamiltonian, 
^o, was used as a zero-order operator. This leads to the perturbation series of 
the type first discussed through second-order by M0ller and Plesset. 55 ' 56 However, 
it is clear that any operator X obeying the relation 

[J?o, = => * = s |*X*I ?|*X*| (55) 

X 

where |«> is an eigenfunction of J^ , may be used to develop a perturbation 
series. The operator 

^shifted = S |kX«|J^I*X*I (56) 

X 

gives rise to the shifted, or Epstein-Nesbet, 57-59 perturbation series. The resulting 
perturbation expansion has the same diagrammatic representation as that 
discussed above for that based on the Hartree-Fock model hamiltonian. The 
corresponding algebraic expressions are as given above in equations (50)-(54) 
except that 

(i) in third-order terms we omit the diagonal scattering terms, i.e. the sum- 
mations are modified as follows 

S S — ► 2 S 

y abed ij abed 

ab^cd 

s s — > z s 

ijk abc ijk abc 

i=£k a^c 

s s — ► s s 

ij¥=kl 

(ii) the denominators are 'shifted' by 

rflJaS = + f abab + J iaai + <fibbi + J^jaa) + J^jbbj (58) 

The use of shifted denominators may also be interpreted as the inclusion of 
certain higher-order terms in the perturbation series based on the Hartree-Fock 
model hamiltonian. 

The sum of the perturbation expansion to infinite order is, of course, in- 
dependent of the choice of the zero-order operator. Assuming that the perturba- 
tion series have converged, the model, or M0ller-Plesset, and shifted, or Epstein- 
Nesbet, perturbation series will give identical results at sufficiently high order. 

Generalizations. — The expansion of Brueckner 10 and Goldstone 11 can only be 
applied when the unperturbed wave function can be described by a single Slater 

55 Chr. Moller and M. S. Plesset, Phys. Rev., 1934, 46, 618. 

66 J. S. Binkley and J. A. Pople, Int. J. Quantum Chem., 1975, 9, 229. 

57 P. S. Epstein, Phys. Rev., 1926, 28, 695. 

58 R. K. Nesbet, Proc. R. Soc. London, Ser. A, 1955, 230, 312, 322. 

59 P. Claverie, S. Diner, and J. Malrieu, Int. J. Quantum Chem., 1967, 1, 751. 



(57a) 
(57b) 
(57c) 
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determinant. Brand ow 51,60 has derived a quasi-degenerate many-body perturba- 
tion theory which has been applied by Kaldor and his co-workers 61 - 62 to some 
open-shell systems which cannot be described by a single determinant in zero- 
order. The formulation is developed with respect to a multi-determinant reference 
function. In Brandow's formalism, 51 ' 60 the occupied orbitals are divided into 
'core' orbitals, which are occupied in all of the determinants included in the 
reference function, and 'valence' orbitals, which are only occupied in some of 
the determinants in the reference function. The problem of applying the many- 
body perturbation theory to open-shell and multiconfiguration reference func- 
tions has been considered by a number of authors. 51 ' 60-74 Brandow's deriva- 
tion 61 ' 60 is time-independent. The derivation of Kuo et al. n is a generalization 
of the time-dependent treatment of Goldstone. 11 Johnson and Baranger 70 give a 
detailed discussion of the so-called folded diagrams which arise in the generalized 
expansion. The work of Lindgren, 72 Levy, 73 and Kvasnicka 74 in this area should 
also be noted. The work of Hegarty and Robb 75 using quasi-degenerate Rayleigh- 
Schrddinger perturbation theory is also of interest. In almost all derivations of 
quasi-degenerate many-body perturbation theory given to date, it is assumed that 
the model space, the space spanned by the determinants in the reference function, 
is complete; that is, all possible occupations of the valence orbitals are included. 
However, this definition of the model space is likely to cause the appearance of 
intruder states. Intruder states are determinants which have to be included in the 
reference function in order to complete the model space but which have an energy 
significantly above that of other determinants in the reference function. The 
presence of intruder states can impair or even destroy the convergence of the 
perturbation expansion. A recent theoretical development in this area is the work 
of Hose and Kaldor 76 which allows an incomplete model space to be used and 
thus enables intruder states to be omitted from the reference function. 

A second generalization of the many-body perturbation theory, which would 
certainly be useful in the calculation of potential energy curves and surfaces, 
would be to use a reference function constructed from non-orthogonal orbitals. 



80 B. H. Brandow, Adv. Quantum Chem., 1977, 10, 187. 
61 U. Kaldor, J. Chem. Phys., 1975, 62, 4364; 1975, 63, 2199. 
" P. S. Stem and U. Kaldor, /. Chem. Phys., 1976, 64, 2002. 
83 C. Bloch, Nuclear Phys., 1958, 6, 329. 

64 J. Des Cloizeaux, Nuclear Phys., 1960, 20, 321. 

65 H. Primas, Helv. Phys. Acta, 1961, 34, 1961. 

66 H. Primas, Rev. Mod. Phys., 1963, 35, 710. 

67 T. Morita, Progr. Theor. Phys., 1963, 29, 351. 

68 H. P. Kelly, Phys. Rev., 1966, 144, 39. 

88 P. G. H. Sandars, Adv. Chem. Phys., 1969, 14, 365. 

70 M. B. Johnson and M. Baranger, Ann. Phys. (N.Y.), 1971, 62, 172. 

71 T. T. S. Kuo, S. Y. Lee, and K. F. Ratcliffe, Nuclear Phys., 1971, A176, 65. 

72 I. Lindgren, /. Phys. B: Atom. Mol. Phys., 1974, 7, 2441; Int. J. Quantum Chem. Symp., 
1978, 12, 33. 

73 B. Levy, in 'Proceedings of the Fourth Seminar on Computational Methods in Quantum 
Chemistry, 1978', ed. B. Roos and G. H. F. Diercksen, Max-Planck-Institut fur Physik und 
Astrophysik, Munich. 

74 V. Kvasnicka, Czech. J. Phys., 1974, B24, 605; Adv. Chem. Phys., 1977, 36, 345. 

75 D. Hegarty and M. A. Robb, Mol. Phys., 1979, 37, 1455. 

76 G. Hose and U. Kaldor, /. Phys. B: Atom. Mol. Phys., 1979, 12, 3827. 
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The work of Gerratt, 77 Gallup, 78 and Goddard 79 suggests a possible reference 
function for use in such a scheme, e.g. ref. 80. Newman 81 ' 82 has examined a 
many-body perturbative formalism which uses biorthogonal sets of orbitals. The 
use of non-orthogonal orbitals has also been discussed by Kvasnicka, 83 Moshin- 
sky and Seligman, 84 Gouyet, 85 Cantu et a/., 86 and Kirtman and Cole. 87 

3 The Algebraic Approximation 

General Remarks. — The algebraic approximiation is fundamental to most 
applications of the methods of quantum mechanics to molecules. When employ- 
ing non-variational techniques, such as perturbation theory, in the study of 
correlation energies it is important to be very clear about the origin of truncation 
effects in a particular calculation. In molecular studies using diagrammatic 
perturbation theory there are, of course, two truncation effects, that arising from 
basis set restrictions (i.e. by invoking the algebraic approximation), and that 
arising from truncation of the perturbation series. In applications of perturbation 
theory to molecules, it is important that one does not obtain good results by a 
fortuitous cancellation of these two sources of error. 

In the first part of this section, the relationship between the solution of the 
Schrodinger equation and the hamiltonian in the space generated by a given basis 
set is discussed in some detail. Since basis set limitations appear to be one of the 
largest sources of error in most present day molecular calculations, the concept of 
a universal even-tempered basis set is discussed in the second part of this section. 
This concept represents an attempt to overcome the incomplete basis set problem, 
at least for diatomic molecules. Further aspects of the basis set truncation prob- 
lem are discussed in the final part of this section. 

The Algebraic Approximation. — The determination of the electronic structure of 
atoms and molecules containing N electrons involves the evaluation of an 
appropriate eigenvalue and associated eigenfunction of a semi-bounded self- 
adjoint hamiltonian operator, J?, in Hilbert space A. A tractable scheme for 
solving such equations is the algebraic approximation in which eigenf unctions 
are parameterized by expansion in a finite set of functions. Integro-differential 
Hartree-Fock equations thus become algebraic equations for the expansion coef- 
ficents. The difficulties of the atomic correlation problem are compounded in 
molecular studies by the impossibility of factorizing the variables in a molecular 
field and use of the algebraic approximation becomes almost obligatory. 

77 J. Gerratt, Adv. Atom. Mol. Phys., 1971, 7, 141 ; and in 'Theoretical Chemistry', ed. R. N. 
Dixon (Specialist Periodical Reports), The Chemical Society, London, 1974, Vol. 1, p. 60. 

78 R. C. Morrison and G. A. Gallup, J. Chem. Phys., 1969, 50, 1214. 

79 W. A. Goddard, Phys. Rev., 1967, 157, 81. 

80 S. Wilson and J. Gerratt, Mol. Phys., 1975, 30, 777; S. Wilson, Mol. Phys., 1978, 35, 1. 

81 D. J. Newman, /. Phys. Chem. Solids, 1969, 30, 1709. 

82 D. J. Newman, /. Phys. Chem. Solids, 1970, 31, 1143. 

83 V. Kvasnicka, Chem. Phys. Lett., 1977, 51, 165. 

84 M. Moshinsky and T. H. Seligman, Ann. Phys. (N.Y.), 1971, 66, 311. 

85 J. F. Gouyet, /. Chem. Phys., 1973, 59, 4637. 

88 A. A. Cantu, D. J. Klein, F. A. Matsen, and T. H. Seligman, Theor. Chim. Acta, 1975,38, 
341. 

« 7 B. Kirtman and S. Cole, /. Chem. Phys., 1978, 69, 5055. 
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The algebraic approximation results in the restriction of the domain of the 
operator to a finite dimensional subspace, Sf, of the Hilbert space d. The 
algebraic approximation may be implemented by defining a suitable ortho- 
normal basis set of M (<N) one-electron spin orbitals and constructing all 
unique iV-electron determinants |/*> using the M one-electron functions. The 

number of unique determinants that can be formed is r\ - ^yj> an< * *l iS tne 

dimension of the subspace £f spanned by the set of determinants. The algebraic 
approximation restricts the domain of & to this ^-dimensional subspace. 

Within the algebraic approximation, the Schrodinger equation may, in prin- 
ciple, be solved by the method of configuration mixing. The wave function is 
then expressed as a superposition of configurations | 11) with linear coefficients 
C A , the optimal expansion coefficients determined by the variation principle. In 
practice, difficulties arise in setting up and solving secular equations of high 
order. Thus, only a small subset of the | is usually employed in the expansion. 

Diagrammatic many-body perturbation theory may also be formulated within 
the algebraic approximation. The domain of the perturbation theory operators 
is restricted to the ^-dimensional space spanned by the | fi) and consequently 
the perturbation theory wave function is generated in terms of an ^-dimensional 
representation. The results of many-body perturbation theory calculations, when 
carried through infinite order, are identical to those of configuration mixing if 
the same basis set, i.e. subspace £f, is used in both calculations. However, in 
practice, the perturbation series must be truncated at some finite order and in 
this case the degree of agreement with the configuration interaction result is a 
measure of the convergence of the perturbation series at this order. 

Let us briefly discuss the relationship between approaches which use basis 
sets and thus have a discrete single-particle spectrum and those which employ 
the Hartree-Fock hamiltonian, which has a continuous spectrum, directly. 
Consider an atom enclosed in a box of radius R, much greater than the atomic 
dimension. This replaces the continuous spectrum by a set of closely spaced 
discrete levels. The relationship between the matrix Hartree-Fock problem, which 
arises when basis sets of discrete functions are utilized, and the Hartree-Fock 
problem can be seen by letting the dimensions of the box increase to infinity. 
Calculations which use discrete basis sets are thus capable, in principle, of 
yielding exact expectation values of the hamiltonian and other operators. In 
using a discrete basis set, we replace integrals over the continuum which arise in 
the evaluation of expectation values by summations. The use of a discrete basis set 
may thus be regarded as a quadrature scheme. 

Once the algebraic approximation has been invoked there is essentially no 
difference between the atomic problem and the molecular problem, except that 
the multicentre integrals which arise in the latter case are more difficult to evaluate. 

Universal Basis Sets. — Historically, it has been necessary to restrict the size of 
basis sets employed in molecular calculations to a reasonably small number of 
functions in order to keep the computation tractable. However, to achieve high 
accuracy, moderately large basis sets are ultimately required, especially if a 
significant fraction of the molecular correlation energy is to be recovered. Since 
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the flexibility of a basis set generally increases as the number of functions is 
increased, the need to optimize parameters becomes less important. This has led 
to the suggestion that a single moderately large basis set, with a reasonable 
choice of parameters, might be suitable for a variety of systems. Such a basis set 
has been termed a 'universal basis set'. 88-92 

It has been shown that Slater type orbitals restricted to Is, 2p, 3d, . . . functions 
can lead to good results for total energies. Orbital exponents, £, may be chosen 
by the 'even-tempered' formula 93-80 

£ fc = k = 1,2, ...,M (59) 

It can be shown that the metric matrix, S, then has the property 

Si,j = Si+i, j+i (60) 

The parameters oc, p, and M for the universal even-tempered basis set thus 
generated have been chosen quite arbitrarily using the following guidelines: 

(i) a must be small enough to ensure a wide range of orbital exponent values, 

(ii) /S must be large enough to avoid possible near linear dependence, 

(iii) M must be large enough to generate a 'near complete' set. 

The following values of a, /?, and M have been shown to be useful in studies 
of atoms and diatomic molecules: 88-92 

Is: oc = 0.5 ft = 1.55 M — 9 

2p: oc = 1.0 p = 1.60 M = 6 (61) 
3d: <x = 1.5 p = 1.65 M = 3 

Thus the orbital exponents of the Is functions range from 0.5 to 16.7 bohr -1 , the 
2p exponents range from 1.0 to 10.5 bohr -1 , and the 3d exponents from 1.5 to 
4.1 bohr -1 . The energy is generally found to be more sensitive to the particular 
choice of P than the choice of a. 

Universal even-tempered basis sets have been employed in calculations on 
first-row and second-row atoms using the matrix Hartree-Fock model. 88 ' 89 
Atomic calculations including electron correlation have also been reported. 90 
For molecules calculations using universal even-tempered basis sets of s, p, and d 
functions have been reported at both the Hartree-Fock level 91 and including 
correlation. 92 For the nitrogen, carbon monoxide, and boron fluoride molecules 
about 80 % of the empirical correlation energy is recovered by taking the perturba- 
tion series through third-order. 

88 D. M. Silver, S. Wilson, and W. C. Nieuwpoort, Int. J. Quantum Chem., 1978, 14, 635. 

89 D. M. Silver and W. C. Nieuwpoort, Chem. Phys. Lett., 1978, 57, 421. 

90 D. M. Silver and S. Wilson, /. Chem. Phys., 1978, 69, 3787. 

91 S. Wilson and D. M. Silver, Chem. Phys. Lett., 1979, 63, 367. 

" S. Wilson and D. M. Silver, J. Chem. Phys., 1980, 72, 2159; unpublished work. 
98 C. M. Reeves and J. Harrison, /. Chem. Phys., 1963, 39, 11. 

94 K. Ruedenberg, R. C. Raffenetti, and R. D. Bardo, in 'Energy Structure and Reactivity', 
Proceedings of the 1972 Boulder Conference on Theoretical Chemistry, Wiley, New York 
1973. 

95 R. C. Raffenetti, /. Chem. Phys., 1973, 59, 5936. 

96 R. D. Bardo and K. Ruedenberg, /. Chem. Phys., 1973, 59, 5956, 5966. 

97 R. C. Raffenetti and K. Ruedenberg, /. Chem. Phys., 1973, 59, 5978. 

98 R. D. Bardo and K. Ruedenberg, /. Chem. Phys., 191 A, 60, 918. 

99 R. C. Raffenetti, Int. J. Quantum Chem. Symp., 1975, 9, 289. 
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Let us comment at this point on the use of the even-tempered method for 
generating orbital exponents. A universal basis set clearly need not be an even- 
tempered one. However, even-tempered basis functions do span the one-electron 
space in a fairly uniform fashion. Unlike a set of exponents determined by 
energy optimization for a particular system, we expect that the even-tempered 
functions will be more suitable for transferring between systems. Even-tempered 
orbital exponents are defined through the recursion 

Cfc+i = t*0 (62) 

A possibly useful generalization of this formula is the recursion 

Cfc+i = C*(/J + ky) (63) 

where y is an additional parameter. 

Several advantages would accrue to the use of a universal basis set. Since most 
electronic structure studies of molecules begin with the valuation of one-electron 
and two-electron integrals over the primitive basis functions, then, for a given set 
of nuclear positions, these integrals could be computed once and used for all 
subsequent studies without regard to the identity of the constituent atoms. This 
transferability extends, of course, to all integrals including multicentre two- 
electron integrals. In order to be an acceptable concept, the universal basis set 
must be capable of providing a uniform description of a series of atoms. Since 
such a basis set is necessarily moderately large, it should yield reasonably high 
accuracy as well as uniformity. 

The concept of a universal gaussian basis set should be briefly discussed since 
gaussian functions are widely used in molecular studies, especially in studies of 
polyatomic systems. Clearly, gaussian basis sets could be transferred from 
molecule to molecule in precisely the manner we have demonstrated for Slater 
functions. However, it is generally accepted that Slater basis functions provide a 
more realistic description of the orbitals especially near nuclei. Integrals over 
universal basis sets are generated once and stored for future use and it would 
therefore seem more profitable to use Slater functions which would generate a 
smaller integral list, especially if a significant fraction of the correlation energy is 
to be recovered. Raffenetti 99 has investigated the use of even-tempered gaussian 
basis sets. His results demonstrate the inferior nature of gaussian functions in 
accurate work. However, it should be noted that multicentre integrals over 
gaussian functions can be evaluated with greater precision than integrals over 
Slater functions. 

It should be noted that the concept of a universal basis set could prove useful 
in almost all molecular studies in which finite basis sets are employed, e.g. 
electron-molecule scattering. Such basis sets lead to accurate results because 
they are moderately large. Furthermore, they afford a certain degree of uni- 
formity in that they are not optimized for one particular property. 

Basis Set Truncation. — Basis set truncation does appear to be one of the main 
sources of error in the majority of calculations on small molecules. The efficiency 
of algorithms based on the diagrammatic perturbation expansion will allow 
increasingly large basis sets to be employed in molecular studies during the next 
few years. 
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Calculations for atomic systems have shown that almost all of the correlation 
energy may be recovered in second-order, if there are no near-degeneracy 
effects in the reference spectrum, when a sufficiently large basis set is employed. 
For the ground state of the neon atom, Eggarter and Eggarter 100-103 have shown 
that 98% of the empirical correlation energy can be recovered by including 
functions with / quantum number up to 6 in the basis set. Since there is essentially 
no difference between the atomic and the molecular problems once a finite basis 
set is introduced (excepting the difficulty of evaluating the multicentre integrals, 
of course), it is expected that similar results will be obtained for molecules once 
sufficiently large basis sets are utilized. 

In order to discuss the convergence of atomic and molecular calculations with 
respect to basis set size, one must have an orderly procedure for extending the 
basis set. This is necessary if the relationship between various restricted basis sets 
and complete bases is to be understood. The recent work of Ruedenberg and his 
co-workers 104 ' 106 is of great interest in this repect. They have examined the 
effective convergence of orbital bases through systematic sequences of gaussian 
primitives. They develop various sets of even-tempered gaussian functions-such 
that 

OL = a(M) (64) 

and 

fi = KM) (65) 

where 

<x(M) -y /S(M) -> 1 W(M)] M -> oo (66) 

as 

M -> oo (67) 

Clearly, this approach can also be used in the case of Slater basis sets and, more- 
over, in the case of universal Slater basis sets. Ruedenberg and co-workers 104 * 105 
have shown that, within the molecular orbital approximation, this systematic 
approach gives a series of energy values which smoothly approach the Hartree- 
Fock limit. Similarly smooth convergence is to be expected in the calculation of 
correlated wave functions and expectation values, and will be the subject of 
future studies in this area. 106 

4 Truncation of the Many-body Perturbation Expansion 
General Remarks. — Perturbation theory forms the basis of a unique approach to 
the calculation of accurate expectation values in that it provides a clearly defined 
order parameter indicating the relative importance of various terms. This order 

100 E. Eggarter and T. P. Eggarter, /. Phys. B: Atom. Mol. Phys., 1978, 11, 1157. 

101 E. Eggarter and T. P. Eggarter, /. Phys. B: Atom. Mol. Phys., 1978, 11, 2069. 

102 E. Eggarter and T. P. Eggarter, /. Phys. B: Atom. Mol. Phys., 1978, 11, 2969. 
10S T. P. Eggarter and E. Eggarter, /. Phys. B: Atom. Mol. Phys., 1978, 11, 3635. 

104 D. F. Feller and K. Ruedenberg, Theor. Chim. Acta, 1979, 52, 231. 

105 M. Schmidt and K. Ruedenberg, /. Chem. Phys., 1979, 71, 3951. 

108 S. Wilson, Theor. Chim. Acta, 1980, 57, 53; ibid., 1980, 58, 31; ibid., 1981, in the press. 
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parameter provides an unambiguous criterion for truncating expansions for 
expectation values. Brandow 107 thus comments that perturbation theory 'seems 
to be the "least biased" of all the many-body techniques'. 

In this section various aspects of the truncation of the perturbation expansion 
are discussed. The use of Pade approximants in perturbation theory is outlined 
and the special invariance properties of the [N+l/N] Pade approximants in 
Rayleigh-Schrodinger perturbation theory are emphasized. Scaling of the zero- 
order hamiltonian is considered and the construction of upper bounds to the 
energy described. Fourth-order and higher-order terms are presented and dis- 
cussed in some detail. Investigations of the effects of quasi-degeneracy on the 
convergence properties of the perturbation series are reviewed. Finally, since 
most other approaches to the electron correlation problem in molecules may be 
regarded as different ways of truncating the many-body perturbation series, a 
brief description of the relationship between many-body perturbation theory and 
some other methods is given. 

Pade Approximants and Perturbation Expansions. — The [P/Q] Pade approxi- 
mant 108 - 109 to a Taylor series T(x) of order M is defined by 

[P/Q](x) = ^(x)/ <2(x); P + Q = M (68) 

where & is a polynomial of order P and j2 is a polynomial of order Q such that 

&{x) - 2(x) [PJQ] (x) = &(x M+1 ) (69) 

where &(x M+1 ) denotes terms of order M+ 1 and higher. 

Pade approximants provide a useful representation of perturbation series. 
They often converge when the usual power series diverges. Unlike the power 
series, Pade approximants can handle a class of function with various types of 
singularities still providing correct uniform convergence. 109 In the Lennard- 
Jones Brillouin Wigner perturbation theory, Pade approximants can be used to 
find bounds to the exact energy. 110 On the other hand, in Rayleigh-Schrodinger 
perturbation theory, the [N+l/N] Pade approximants are special 111 - 112 in that, 
when the expansion parameter is set equal to unity, the numerical value of this 
approximant is invariant to two modifications in the zero-order hamiltonian 
operator, namely a change of scale and a shift of origin. 

Let us consider these two modifications of a given zero-order operator, J^ ; 
that is a uniform displacement of the zero-order energy spectrum and a uniform 
change of scale in the spacing of the zero-order energy levels. Consider the zero- 
order operator 

= putfo + vt (70) 

where / is the identity operator and n and v are scalars. The corresponding 
perturbation operator is 

j#W = + (l-^o-vl (71) 

107 B. H. Brandow, in 'Effective Interactions and Operators in Nuclei', ed. B. R. Barratt, 
Springer- Verlag, Berlin, 1975. 

108 H. Pade, Ann. sci. Ecol. norm. sup. Paris (Suppl.), 1892, 9, 3. 

108 G. A. Baker, 'Essentials of Pade Approximants', Academic Press, New York, 1975. 

110 O. Goscinski, Int. J. Quantum Chem., 1967, 1, 769. 

111 S. Wilson, D. M. Silver, and R. A. Farrell, Proc. R. Soc. London, Ser. A, 1977, 356, 363. 
118 E. Feenberg, Ann. Phys. (N. Y.\ 1958, 3, 292. 
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Note that when the perturbation parameter is put equal to unity the total 
hamiltonian operator is recovered, i.e. 

3# = (j*^,* + Aj^.»0a-i (72) 

The energy coefficients, Ef> v , where i denotes the order of perturbation, may be 
related to the values of E\>° by the expressions 

E£ v = fiE\ -o + v (73) 

E* v = El'° + (l-fi)E{'°-v (74) 



and 



E "n V = ^% (j-z) <>-iy*-qr n = 2,3,... (75) 



It can be shown 111 that the [N+l/N] Pade approximants alone among all 
Pad6 approximants of order 2N+ 1 are invariant to an arbitrary choice of /* 
and v, when A is set to unity. 

The [N+ 1/N] Pad6 approximants to the energy series also have the advantage 
that they have a linear dependence on the perturbation parameter, A, as it 
becomes large. Furthermore, the Goldstone level shift formula 11 may be written 

AE = <0 O \3?i 0\ #o>l (76) 

where the evolution operator, C, is unitary if all terms are included but not if 
the expansion is truncated. A unitary approximation to can be obtained by 
forming the [N/N] operator Pad6 approximant. 113 This has the same form as a 
function of the perturbation parameter as the |W+ 1/iV] scalar Pade approximant 
to the level shift. 

From third-order perturbation calculations, the [2/1] Pade approximants can 
be formed. Several additional reasons can be advanced for their use: 

(i) there is a large amount of evidence demonstrating that this leads to im- 
proved results, 114 

(ii) this may be regarded as the choice of scaling parameter which makes the 
third-order energy vanish, 115-117 

(iii) it may also be regarded as the energy expression for which the energy is a 
minimum with respect to choice of scaling parameter. 118 

It should be noted that for one-electron properties, such as dipole moments, 
the [N/N] Pade approximants are invariant to changes of scale and shifts of 
origin in the reference spectrum 8 whereas for second-order properties, such as 
polarizabilities, the [N/N+ 1] Pad6 approximants are to be preferred. Indeed, 
for polarizabilities the use of the form 

has been advocated by a number of authors. 8 ' 119 - 120 

113 J. L. Gammel and F. A. McDonald, Phys. Rev., 1966, 142, 1245. 

114 S. Wilson, /. Phys. B: Atom. Mol. Phys., 1979, 12, 1623. 

115 A. T. Amos, J. Chem. Phys., 1970, 52, 603. 

A. T. Amos, Int. J. Quantum Chem., 1972, 6, 125. 
ii? A. T. Amos, /. Phys. B: Atom. Mol. Phys., 1978, 11, 2053. 
H8 d. T. Tuan, Chem. Phys. Lett., 1970, 7, 115. 
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Pad6 approximates to perturbation series have been considered by numerous 
authors, 121 particularly those interested in effective interactions in nuclei. They 
have been used in studies of the 'intruder state' problem which can arise when a 
multideterminantal function is used as a reference in a perturbative scheme. 122 

Scaling of the Zero-order Hamiltonian. — It has been demonstrated 123-128 that 
improved results can be obtained from second-order Rayleigh-Schrodinger 
perturbation studies of electron correlation energies using very large basis sets 
by introducing a scaling parameter, in the zero-order hamiltonian, 3& . The 
scaling parameter may be determined, from calculations taken to third-order in 
the energy and using a basis set of moderate size, by imposing the condition that 
the third-order energy coefficient corresponding to the zero-order operator, 
fi-^o, be zero. The correlation energy is then approximated by iir x E^ where E 2 
is the second-order energy resulting from the calculation using the large basis set. 

Second-order calculations with large basis sets are computationally tractable. 
One of the most time consuming parts of many calculations is the transformation 
of the two-electron integrals over atomic functions to integrals over molecular 
orbitals. For second-order calculations only a restricted list of two-electron 
integrals is required (although the restricted transformation is still an n 5 process, 
where n is the number of basis functions). It is envisaged that the concept of a 
universal even-tempered basis set will prove useful in developing the large basis 
sets required for these calculations. 88-92 The scaling of the zero-order hamiltonian 
is also useful in the calculation of molecular properties and has been demon- 
strated 120 for the polarizability of the hydrogen molecule. 128 ' 127 

Modified Potentials. — It is possible to modify the perturbation series by using 
alternative potentials in the zero-order hamiltonian. The modified potential 

W-H'i) = N i 1 f dr 2 <f>*( n )r$ (1 -^12) M'*> (78) 

has been widely used in atomic studies. 6 Recent calculations 128 ' 129 on the 
hydrogen fluoride molecule have indicated that there is little difference between 
results obtained using the Hartree-Fock potential and the above modified 
potential when all terms are included through third-order. For the FH molecule 
using unshifted denominators, 80.2% of the empirical correlation energy was 

119 G. Howat, M. Trsic, and O. Goscinski, Int. J. Quantum Chem., 1977, 11, 283. 
"OS. Wilson, Mol. Phys., 1980, 39, 525; S. Wilson and A. J. Sadley, to be published. 
121 G. A. Baker, Adv. Theor. Phys., 1965, 1, 1 ; O. Goscinski and E. Brandas, Phys. Rev. A, 1970, 
1, 552. 

181 E. M. Krenciglowa and T. T. S. Kuo, Nuclear Phys., 1974, A235, 171 ; H. M. Hoffmann, 
S. Y. Lee, J. Richert, and H. A. Weidenmuller, Ann. Phys., 1974, 85, 410; J. M. Leinaas 
and T. T. S. Kuo, Ann. Phys., 1978, 111, 19; T. H. Schucan and H. A. Weidenmuller, 
Ann. Phys., 1972, 73, 108; J. Richert, T. H. Schucan, M. H. Scrubel, and H. A. Weiden- 
muller, Ann. Phys., 1976, 96, 139; T. H. Schucan and H. A. Weidenmuller, Ann. Phys., 
1973, 76, 483. 

188 S. Wilson, J. Phys. B: Atom. Mol. Phys., 1979, 12, L135. 

124 F. W. Byron and C. J. Joachain, /. Phys. B: Atom. Mol. Phys., 1979, 12, L597. 

128 S. Wilson, /. Phys. B: Atom. Mol. Phys., 1979, 12, L599. 

"«T. Itagaki and A. Saika, Chem. Phys. Lett., 1977, 52, 530. 

1.7 T. Itagaki and A. Saika, /. Chem. Phys., 1979, 70, 2378. 

1.8 D. M. Silver and R. J. Bartlett, Phys. Rev. A, 1976, 13, 1. 

1.9 D. M. Silver, S. Wilson, and R. J. Bartlett, Phys. Rev. A, 1977, 16, 477. 
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recovered in second-order when the Hartree-Fock potential was employed, 
whereas 101.2 % of the correlation energy was recovered when the V 1 *- 1 modified 
potential was used. The [2/1] Pade approximant to the perturbation series using 
the Hartree-Fock potential gave 79.6 % of the correlation, which is close to the 
77.8% given by the [2/1] Pad6 approximant to the series based on the modified 
potential. The infinite-order result should, of course, be independent of the choice 
of zero-order hamiltonian. 

Upper Bounds to Total Energies. — In performing a calculation of the energy to 
third-order one also implicitly calculates the wave function to first-order. Let the 
wave function through first-order be written as 

#o + Mi (79) 

where the parameter A has been introduced. The truncated expansion for the 
energy is not bounded : however by substituting the first-order wave function in 
the Rayleigh quotient, an upper bound is obtained. This may be written as 32 

£ var (A) = Eo + Ei + (n ~P + %Zf E3 (80) 

where the overlap integral, A n , is given by 

Au = iss(^) 2 (81) 

A may be regarded as a variational parameter. Its optimal value may be obtained 
by invoking the variation principle, giving 



*o P t = An 1 [VC 2 + An-C] (82) 

where 

f = Hl-W (83) 

Unfortunately, by forming an upper bound one introduces terms depending on 
the square of the number of electrons and thus the quality of this upper bound 
deteriorates as the number of electrons is increased. This is illustrated in Table 1 
for an array of well-separated nitrogen molecules. 



Table 1 Size-consistency error in an array of well separated N 2 Molecules^ 



m 


n a 




e e 


1 


14 


0.00 


-0.45 


2 


28 


0.09 (10%) 


-0.90 


3 


42 


0.30(22%) 


—1.35 


4 


56 


0.59 (33%) 


—1.80 


5 


70 


1.00(44%) 


-2.25 



tBased on results given in ref. 92. 

° Number of electrons. b Fourth-order size-consistency error; percentage of estimated cor- 
relation energy given in parentheses. c Estimated correlation energy. 

Fourth-order and Higher-order Terms. — Most methods currently employed in the 
study of electron correlation in atoms and molecules, for example the coupled- 
electron pair theory, limited configuration mixing, etc., may be regarded as third- 
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order theories in that they neglect or approximate fourth-order and higher-order 
terms in the perturbation analysis of the energy. Diagrammatic many-body 
perturbation theory offers a systematic scheme for extending such calculations. 
The perturbation series forms the basis of a balanced treatment of the correlation 
problem in that it provides a clearly defined order parameter which indicates the 
relative importance of various terms. There is considerable interest in the evalua- 
tion of the fourth-order terms in the many-body perturbative expansion since 
these terms represent, at least in part, the dominant corrections to most tech- 
niques currently being employed in calculating atomic and molecular correlation 
energies. 

Unlike the second-order and third-order energy diagrams, the fourth-order 
diagrams can involve intermediate states which are singly-excited, doubly- 
excited, triply-excited, and quadruply-excited with respect to the Hartree-Fock 
reference function. 8 - 130 

There are four diagrams, using the Brandow convention 51 which involve an 
intermediate state which is singly-excited with respect to the Hartree-Fock 
function, these are shown in Figure 5. Diagrams As and Ds are related by time 
reversal, whereas Bs and Cs are related by complex conjugation. 130 

3 f\~~3t /%3 

3 V-3 V:3 V:3 

A s B s C s D s 

Figure 5 Fourth-order diagrams which involve singly-excited intermediate states 

The twelve diagrams of the fourth-order which involve only doubly-excited 
intermediate states are displayed in Figure 6. Diagrams (Bd, Cd), (Ed, Fd), and 
(Gd, Hd) are related by complex conjugation. Diagrams (Ad, Dd), (Ed, Hd), 
(Fd, Gd), and (Jd, Kd) are related by time reversal. 

The diagrams which involve a triply-excited intermediate state in fourth-order 
are shown in Figure 7. Explicitly, diagrams At, Bt, Or, and Dt, for example, 
correspond to the expressions 

E 4 (At) = -i £( J iiabj akcdf cbekj edi})l(@ijab@ ' ijkbcd^ijde) (84) 

£4(J3t) = -i 2(/ i)ab faked f edejf ebik)l(@Vab3>i)kbcd@ikbe) (85) 

£4(Ct) = -i S (fi)abfklicfmbklfacm))l(3 ! i}ab@ l Wabc@ l )mac) (86) 

£i(Dt)= —i Z(f i} abJ'klicf'mckjf~abmi)l(3'ijab@jklabc&emab) (87) 

Time-reversal symmetry can be seen to relate diagrams At, Bt, Et, and Ft to 
diagrams Dt, Ct, Ht, and Gt, respectively. Diagrams (It, Lt) and (Jt, Kt) form 
complex conjugate pairs. Diagrams (Mt, Pt) and (Nt, Ot) are related by time 
reversal and (Mt, Nt) and (Ot, Pt) are complex conjugates. 

130 S. Wilson and D. M. Silver, Int. J. Quantum Chem., 1979, 15, 683. 
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Id Jd Kd Lq 

Figure 6 Fourth-order diagrams which involve doubly-excited intermediate states 



Finally, there are seven fourth-order terms which involve an intermediate 
state which is quadruply-excited with respect to the Hartree-Fock reference 
function. These diagrams are shown in Figure 8. Time reversal relates diagrams 
(Bq, Cq), (Dq, Gq), and (Eq, Fq). Explicitly, diagrams Aq, Bq, and Cq, for 
example, correspond to the expressions 

2?4(Aq) = £ £ (J'ijabJ'klc<I</cbuJ'adk))l(@t}ab@i]klabcd3kjad) (88) 
ijkl abed 

£4(Bq) = -fc £ £ (J'ijab,fklc<l c fcdljtfabkl)l(@ l i]ab@lWabcd@klab) (89) 
ijkl abed 

Ei(Cq) = £ £ (J'ijabJ'klcdJ'abklJ'cdii)l(3>iJab@iJklabca@i}cti) (90) 
ijkl abed 

The fourth-order quadruple-excitation energy diagrams arise from the dis- 
connected wave function diagrams shown in Figure 9. This observation leads to 
considerable simplification in the evaluation of this energy component. For 
example, using the identity 

(@i)al>3'wiabcd2)ijca)- X + (®i)ab@iWabcd@klab)- 1 = (@ijab@i}cd@klab)- 1 (91) 

the sum of the energy components corresponding to diagrams Bq and Cq may 
be written as 

JE 4 (Bq + Cq) = -A" S (J'ijabJ'klcdJ'cdi)J'abki)K3>ijab@ijcd3'kl a b) (92) 

It should be noted that this simplification does not occur when the shifted, or 
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Figure 9 Connected and disconnected wave function diagrams: (a) second-order connected 
triple-excitation diagram, (b) second-order disconnected quadruple-excitation 
diagram, (c) third-order disconnected triple-excitation diagram, (d) third-order 
connected quadruple-excitation diagram 

Epstein-Nesbet, perturbation scheme is used. All of the quadruple-excitation 
terms may be written in the form 

M% (93) 
f 

where 
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and 

gn — 2a~^ y?->) 



M is a constant. ^ and v denote compound indices representing a maximum of 
four actual indices. The intermediates and g^ differ only in that in the latter 
there is an additional denominator factor. A number of calculations of the 
quadruple-excitation component of the fourth-order energy have been re- 
ported. 131 " 140 

The triple-excitation fourth-order energy, in contrast to the quadruple- 
excitation component, arises from connected wave function diagrams. The 
algorithm required to evaluate this energy component is considerably less 
tractable than that for the quadruple-excitation energy, depending on « 7 , where 
n is the number of basis functions. The triple-excitation diagrams can be written 
in terms of the intermediates. 



and 



fm; ab c = 2 l^-l^ (96) 



_ _ £ jlba^iklc SQ-i\ 

gim«>c = S -3—- (97) 



Diagrams At, Bt, Ct, and Dt correspond to the expressions 

£ 4 (At) = — \ S Hfijkiabcfmiacbl^ijkabc (98) 

£4(Bt) = -i S HfM;abcfik);abcl@>Mabc (99) 

£4(Ct) = —$Y,Xgijk;abcgW,acbl@iJkabc (100) 

E 4 (Dt) = —\ S l£gi)k;abcgik);abcl@'ijkabc (101) 

The computation of these energy components has been rendered reasonably 
tractable by 

(i) using a spin-free formalism; 

(ii) recognizing certain permutational symmetry properties of the inter- 
mediates fijk ;ab c and gi)k ;aba ', 



131 S. Wilson and D. M. Silver, in 'Proceedings of the Fourth Seminar on Computational 
Problems in Quantum Chemistry', ed. B. Roos and G. H. F. Diercksen, Max-Planck- 
Institute fur Physik und Astrophysik, Munich, 1978. 

182 S. Wilson and D. M. Silver, Mol. Phys., 1978, 36, 1539. 

138 S. Wilson and D. M. Silver, Comput. Phys. Comm., 1979, 17, 47. 

134 R. Krishnan and J. A. Pople, Int. J. Quantum Chem., 1978, 14, 91; R. Krishnan, 
M. J. Frisch, and J. A. Pople, /. Chem. Phys., 1980, 72, 4244; M. J. Frisch, R. Krishnan, 
and J. A. Pople, Chem. Phys. Lett., 1980, 75, 66. 

135 L. T. Redmon, G. D. Purvis, and R. J. Bartlett, /. Chem. Phys., 1978, 69, 5386. 

138 R. J. Bartlett, I. Shavitt, and G. D. Purvis, /. Chem. Phys., 1979, 71, 281. 

187 I. Hubad, M. Urban, and V. Kello, Chem. Phys. Lett., 1979, 62, 584; M. Urban, I. Huba£, 
V. Kello, and J. Noga, /. Chem. Phys., 1980, 72, 3378. 

188 S. Wilson and V. R. Saunders, J. Phys. B: Atom. Mol. Phys., 1979, 12, L403 ; corrigenda, 
1980, 13, 1505; M. F. Guest, and S. Wilson, Chem. Phys. Lett., 1980, 72, 49; S. Wilson 
and M. F. Guest, ibid., 1980, 73, 607. 

139 V. Kvasnicka, V. Laurinc, and S. Biskupic, Chem. Phys. Lett., 1979, 67, 81. 

140 D. M. Silver, S. Wilson, and C. F. Bunge, Phys. Rev., 1979, A19, 1375. 
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(iii) using a vector-processing computer. 
Again a number of calculations of the triple-excitation component of the correla- 
tion energy have been reported. 138 - 13B > 141 

Disconnected triple-excitation wave function diagrams and connected quad- 
ruple-excitation diagrams first contribute to the energy in fifth-order. A full set 
of fifth-order diagrams, using the notation of Paldus and Wong 1 "- 1 " for the 
sake of brevity, is given in Figure 10. 9 Calculations through fifth-order would 
enable [3/2] Pade approximants to the energy series to be constructed and 
variational upper bounds calculated from second-order wave functions. Full 
fifth-order calculations do not appear to be possible at present. However, the 




Figure 10 Fifth-order diagrams, using the modified Hugenholtz notation of Paldus and 
Wong 



141 S. Wilson and V. R. Saunders, Comput. Phys. Comm., 1980, 19, 293. 

142 J. Paldus and H. C. Wong, Comput. Phys. Comm., 1973, 6, 1. 

143 H. C. Wong and J. Paldus, Comput. Phys. Comm., 1973, 6, 9. 
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recent work by Kvasnicka et a/. 144 on Wigner's (2«+l) rule in many-body 
perturbation theory will certainly simplify fifth-order studies. 

Quasi-degeneracy Effects. — The convergence properties of the non-degenerate 
formulation of the many-body perturbation theory deteriorate when quasi- 
degeneracy is present in the reference spectrum. In view of its simplicity, however, 
there is considerable interest in exploring the range of applicability of the non- 
degenerate formalism. 

The low-lying 2p state in the Be atom ground state makes this a severe test 
of the non-degenerate perturbation series. It is an appropriate choice of system 
for investigating the limitations of a non-degenerate formulation of perturbation 
theory to systems involving near degeneracy. Full configuration mixing calcula- 
tions and third-order many-body perturbation theory studies of the ground state 
of the beryllium atom have been performed 140 within the same basis sets. For a 
basis set containing only functions of s symmetry, there are no near degeneracy 
effects and the degree of agreement between the configuration mixing and 
perturbative studies is good. For a basis set of functions having s, p, and d 
symmetry, the perturbation series based on the Hartree-Fock model zero-order 
hamiltonian gave only 93 % of the correlation correction given by the configura- 
tion mixing calculation after forming the [2/1] Pade approximant to the energy 
series. On the other hand, by forming the [2/1] Pade approximant to the shifted 
perturbation series, 99.5 % of the configuration mixing correlation energy was 
recovered. 

Third-order non-degenerate many-body perturbation theory studies of the 
ground state of the CH + ion have been reported. 145 Calculations were made 
using the same basis set of Slater functions employed by Green et a/. 146 in their 
configuration mixing calculations, which included all single- and double-excita- 
tions from a two-configuration reference function. For this system quasi-degene- 
racy effects become increasingly important as the nuclear separation is increased. 
It was again demonstrated that [2/1] Pade approximants to the shifted perturba- 
tion series are useful when near degeneracy is present. 

One interpretation of the significance of denominator shifts, which are used in 
the shifted perturbation series, is that they provide a method for summing 
certain higher-order diagonal scattering diagrams in the model perturbation 
scheme to infinite order. 32 In this sense, the shifted third-order results include 
all third-order contributions in the model scheme plus a selection of additional 
diagrams summed to infinite order. The quasi-degeneracy problem can be 
overcome in the Be atom and the CH+ ion for internuclear distances < 4 bohr by 
including certain higher-order diagrammatic terms. Alternatively the shifted 
perturbation scheme can be interpreted as an expansion based on an alternative 
zero-order hamiltonian 147 whose spectrum does not contain near degeneracies. 
We prefer this latter interpretation since it represents a more balanced treatment 

144 V. KvasniCka, V. Laurinc, and S. Biskupie Mol. Phys., 1980, 39, 143. 

145 S. Wilson, /. Phys. B: Atom. Mol. Phys., 1979, 12, 1623. 

148 S. Green, P. S. Bagus, B. Liu, A. D. McLean, and M. Yoshmine, Phys. Rev. A, 1972, 5, 
1614. 

147 P. Claverie, S. Diner, and J. Malrieu, Int. J. Chem., 1967, 1, 751. 
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in that all terms are included through third-order and no partial evaluation of 
higher-order terms is attempted. 

In order to eliminate the arbitrary use of the model and shifted perturbation 
series, depending on the degree of degeneracy in the corresponding reference 
spectra, the following zero-order hamiltonian has been considered: 148 

^0 = (1 - yW).#V odel + i M^g hifted (102) 

where fi is an arbitrary scalar. Prototype calculations have been performed for 
the Be atom in which quasi-degeneracy effects arise, and for the Ne atom in 
which no near degeneracy is present. For both Be and Ne, £[3/0] and £[2/1 ] were 
found to exhibit minima with respect to p (although it should be noted that 
£[3/0] and £[2/1] are not bounded). For the Be atom the [2/1] Pad6 approximant 
to the energy series takes its minimum value at ,u~0.9 while for Ne £[2/1] is a 
minimum at /*~0.2. For Be this minimum occurs closer to /« = 1, that is the 
shifted, or Epstein-Nesbet series, than a* = 0, the model or Meller-Plesset series. 
The reverse situation applies in the case of the Ne atom. The £[2/1] energy 
values for both Be and Ne are less sensitive to the value of fi than the usual power 
series, £[3/0]. 

Clearly, given a set of zero-order functions which are eigenfunctions of some 
zero-order operator, «?^ , it is possible to develop a perturbation with respect to 
any operator, A, provided that [3^ ,A] = 0. The infinite-order perturbation 
expansion should of course be independent of the choice of reference hamiltonian 
and thus the degree of agreement between two perturbation series when truncated 
at some finite order gives a qualitative measure of convergence. Strictly, if two 
perturbation expansions agree exactly at some order this does not mean that 
they have converged, but only that they are converging at the same rate. The 
use of a combination of two operators as described above allows a set of zero- 
order operators and the corresponding perturbation series to be investigated. 

Obviously, when severe near degeneracy is present in the reference spectrum, a 
full quasi-degenerate formulation of the perturbation series should be employed. 
However, there is considerable value in applying the much simpler non-degener- 
ate formalism to as wide a range of problems as possible. 

Comparison with Other Methods. — Most methods for performing accurate cal- 
culations of the electronic properties of molecules involve some finite or infinite 
expansion for the wave function and the corresponding expectation values. This 
is the case in, for example, the method of configuration mixing, the many-body 
perturbation theory, cluster expansions, coupled-pair approximation, etc. Each 
method leads to the exact wave function if all terms in the expansion are evalu- 
ated. In practice, the expansions have to be truncated and methods differ only in 
the manner in which this truncation is effected. However, this truncation can 
significantly affect the theoretical properties, and to some extent the computa- 
tional properties, of a method. The mere fact that one method includes more 
terms than another does not mean that it is superior. Which terms one leaves out 
is just as important as the terms one actually evaluates. Thus in the method of 
configuration mixing limited to single- and double-excitation with respect to a 

" 8 S. Wilson, Chem. Phys, Lett., 1979, 66, 255. 
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single determinantal reference function one includes many terms, corresponding 
to unlinked diagrams in the perturbation analysis, which cancel with terms in- 
volving higher-order excitations. 

Perturbation theory provides a clearly defined order parameter in the expansion 
for expectation values giving a 'least-biased' indication of the importance of 
various terms. To quote Brandow: 107 

'The structure of each formalism tends to suggest that certain types of approxi- 
mations are most reasonable or "natural", regardless of the actual quantitative 
characteristics of the physical system. The many-body literature is full of papers 
where people have tried to make the physics fit into their preconceived approxima- 
tion schemes, instead of vice versa. . . . perturbation theory . . . seems to be the 
"least biased" of all the many-body techniques.' 

The relationship between the coupled-electron pair approximation (c.e.p.a.) 
and the many-body perturbation theory has been discussed in detail by Ahl- 
richs. 149 All of the methods denoted by c.e.p.a. (x) (x=0, 1, 2, 3) may be related 
to the summation of certain classes of diagrams in the many-body perturbation 
theory to infinite order. For example, c.e.p.a. (0), which is Cizek's linear approxi- 
mation or Hurley's c.p.a. (0) ansatz 150 is equivalent to the summation of all 
double-excitation linked diagrams in the perturbation series. This is also denoted 
d.e.m.b.p.t. (double excitation many-body perturbation theory) by some 
workers. 151 - 182 

In the coupled pair approximation, 150 all diagrams which can be formed by 
considering products of disconnected double-excitations are summed to infinite 
order. Thus in fourth-order, the coupled pair approximation includes the linked 
double-excitation and linked quadruple-excitation diagrams shown in Figure 6 
and Figure 8, respectively. The coupled pair approximation may be represented 
diagrammatically as follows: 

-(V._.V)= 1+ V...V^(V...V) 2+ - <■» 

The order of perturbation at which various levels of excitation first arise is 
illustrated in Figure 1 1 for three different reference functions. In Figure 11(a), the 
Hartree-Fock orbitals are used to form the reference function, in Figure 11(b) 
the bare-nucleus model is used in zero-order, while in Figure 11(c) Brueckner 
orbitals are used to construct the reference function. 

The method of single- and double-excitation configuration mixing is not 
size-consistent. The method remains size-inconsistent even if a multi-configura- 
tional reference function is used. For a configuration mixing calculation, with 



149 R. Ahlrichs, Comput. Phys. Comm., 1979, 17, 31. 

150 A. C. Hurley, 'Electron Correlation in Small Molecules', Academic Press, 1976. 

151 R. J. Bartlett and I. Shavitt, Int. J. Quantum Chem. Symp., 1977, 11, 165. 

182 M. R. A. Blomberg and P. E. M. Siegbahn, Int. J. Quantum Chem., 1977, 14, 583. 
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respect to a single determinant, a correction to fourth-order can be made by 
means of Davidson's formula 2 ' 153 ' 154 



where £b is the double-excitation configuration mixing energy and C is the 
coefficient of the reference function. The formula 2 • 131 > 155 - 156 



is correct to fifth-order. 

All of the methods discussed above, and some others which we have not 
discussed, represent a partial evaluation of fourth-order and higher-order terms 
in the many-body perturbation series. This must be done with great care since 
often there is a large degree of cancellation between terms in a given order. For 
example, it has been demonstrated that the two-body and many-body compo- 
nents of the energy are of opposite sign and for a system such as the hydrogen 
fluoride molecule the latter can have a magnitude which is ~ 25 % of the two-body 
component. 167 In a fourth-order calculation on water, it was found that there 
was a significant degree of cancellation between the fourth-order double and 
quadruple- excitation components. 156 



General Remarks. — In this section, the computational aspects of the application 
of the many-body perturbation theory to molecular correlation energies are 
discussed. 

Any molecular calculation starts with the evaluation of integrals over the basis 
functions. This is usually, but not necessarily, followed by a self-consistent field 
calculation and a transformation from integrals over atomic basis functions to 
integrals over molecular orbitals. Full details of these particular phases of cal- 
culation are well documented elsewhere and we do not consider them further here. 

Third-order calculations are considered in the next section. This is followed by 
a brief discussion of the computation of higher-order terms and of the evaluation 
of 'bubble' diagrams which are required when molecular properties are calcu- 
lated or when a reference function other than the closed-shell Hartree-Fock 
function is employed. The impact of the new generation of computers, which 
have vector processing capabilities, on many-body perturbative calculations is 
discussed very briefly in the final section. 

Third-order Many-body Perturbative Calculations.— Starting from a list of one- 
electron and two-electron integrals over molecular orbitals, the evaluation of the 
second-order and third-order energies 12-14 consists of two stages. In the first stage 
the integrals are sorted into various types 12 while in the second stage these integral 
lists are used to compute the required energy components. 13 ' 14 

153 E. R. Davidson, in 'The World of Quantum Chemistry' Proceedings of the First Inter- 
national Congress on Quantum Chemistry, ed. R. Daudel and B. Pullman, D. Reidel, 1974. 

154 S. R. Langhoffand E. R. Davidson, Int. J. Quantum Chem., 1974, 8, 61. 

155 P. E. M. Siegbahn, Chem. Phys. Lett., 1978, 55, 386. 

158 S. Wilson and D. M. Silver, Theor. Chim. Acta, 1979, 54, 83. 
157 S. Wilson and D. M. Silver, /. Chem. Phys., 1977, 66, 5400. 
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In the integral sorting phase of the calculation, the following lists of integrals 
are created: 

(1) <IJ\v\KL> 

(2) <AB\v\CD> 

(3) <IA\v\JK> v = ± (106) 

(4) <IA\v\BCy 

(5) <H\v\AB> 

(6) <IA\v\JB> 

where I,J,K, ... are used to denote occupied orbitals and A,B,C, . . . are used to 
denote virtual orbitals. Integral types (3) and (4) are not required when a Hartree- 
Fock reference function is used in calculations through third-order. It is also 
convenient to produce a list of integrals of the form (PQ \v \ PQ) and <PQ | v | QPy 
which are required for the denominator shift factors and, together with the one- 
electron integrals, to evaluate the matrix Hartree-Fock energy. 

To illustrate the calculation of the energies corresponding to the four Brandow 
diagrams which arise through third-order we shall consider the third-order 'hole- 
particle' diagram in some detail. 14 The various spin types which can arise are 
shown in Figure 1 2. The first step is to process the integrals of the type <// 1 v \ A By . 
For an atom or molecule which is described by JV doubly-occupied orbitals, N 
secondary lists are created. The 7th of these secondary lists contains the integrals 
(IJ\v\AB} for all / and A^B, in order. The algorithm is 'integral driven.' 
Integrals of the type </C| v \AK} are read into high-speed core, for all A ^ C and 
a given I^K. Next the integrals <IJ\v\AB} and (JK\v\BCy are read into core 
from the 7th and Kth secondary lists, respectively. All integrals depending on a 
given I, J, and K are in core simultaneously. The integrals of the form <//| v\AB} 
and (JK\v\BCy are read for all values of J in turn until the 7th and ^Tth secon- 
dary lists have been read through completely. A new set of integrals <JC\v\AKy 
is then processed. For a given set of integrals depending on /, /, and K the 
following spin-free intermediates are formed: 

F i IVAn _ y «// \v\ABy-<IJ\v \BA»«JK 1 v \BC> -<JK\v\CB» nfm 

~ | 9IJAB&JKBC (107) 

FhKAC = S W\'\A »^\ *\BAyKJK\v\BC> ( 

f$jkac = a ^ZMABl WvlBc y 

B 3)ijab9Ijkbc 
FW = E mi^M^h^^bd^) (no) 

B ^IJAB^JKBC V 

I^S™ <>») 
The two remaining particle states can then be summed over to give 

S? JK = £ FJ&ac «/C \v\AKy-iIC\v\ KAy) (112) 

AC 
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13 14 
Figure 12 Spin types which arise for the third-order ring diagram 



8?)k = S {FWkac «IC\v \AK>~<IC \v [KA» + Fj? KAC <IC \v \KA>} (113) 

AC 



8?} K = E F^ac^CUIAK} 

AC 

8 { $k= X Ftf KAC <IC\v\AK> 

AC 



(114) 
(115) 

(116) 



Finally, the third-order 'hole-particle' energy is given by 

£ 3 (hp) = - 2^ QtfJ* + 8?)k + 8?}k + 8Wk) 

Similar algorithms have been devised to evaluate the other third-order and the 
second-order energy components. 

Higher-order Terms. — Diagrammatic perturbation theory provides a tractable 
scheme for calculating the dominant components of the correlation energy which 
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may be associated with triple-excitations and quadruple-excitations. As an 
example we consider the evaluation of the quadruple-excitation fourth-order 
linked diagram energy component. 

When the matrix Hartree-Fock operator is used as a reference hamiltonian, 
that is the scheme of M0ller and Plesset, the evaluation of the quadruple- 
excitation energy component is particularly easy. Only one type of integral is 
required, namely the (IJ\v\ABy integrals. For a system which is described by N 
doubly occupied orbitals, N secondary lists are created in the manner described 
in the previous part of this section for the third-order hole-particle energy. The 
computation of the correlation energy corresponding to the diagrams shown in 
Figure 8 requires that certain blocks of integrals be in computer core storage simul- 
taneously. This is achieved by two file handling schemes; one for diagram Aq 
and the sum of diagrams Eq and Gq, and one for the sum of diagrams Bq and 
Cq, and Dq and Eq. In our method, diagram Aq and the sum of Eq and Gq are 
evaluated by employing the same file handling scheme. Blocks (JJ) and (JK) are 
placed in core storage simultaneously. This is done by reading simultaneously 
through the 7th and Kth. secondary lists for all values of /. The 7th and Kih 
secondary lists are read through completely for each value of / and K. When 
computing the energy component corresponding to the sum of diagrams Bq and 
Cq, we require blocks (JJ) and (KL) in core storage at the same time. This is 
achieved by reading the block (//) from the primary list of integrals, rewinding 
and reading all (KL) less than (//). The next value of (//) can then be treated, the 
file rewound and the whole process repeated. 

Spin orthogonalities considerably simplify the computations. For example, in 
the evaluation of diagram Aq six spin-free intermediates can be defined. Three 
of these take the form 

«// \v\AB>-<IJ\v\ BA»«KJ \v \AD>-(KJ\v \DA» 

_v + <IJ\v \BAXKJ\ v\DA> , nTi 

JMbd - L WJJH K } 

JA 

«// \v\AB>-<IJ\v \BAy)<KJ | v\DA> 
, (2J _ v + (IJ\v\BA>((KJ\v\ADy- <KJ\v\DAy) am 
JlK,SD ~ ~ \i Quad U ' 

J A 



_ V <U\v\ABXKJ\v\AD > am 



JA 



while the remaining three merely include an extra denominator factor [see 
equation (95)]. The energy corresponding to diagram Aq is then given by 

£ 4 (A Q ) = S S {flSBDg&BD + WbdS^BD + m BD g®BD} ( 120 ) 
IK BD 

We note that iterative schemes have been devised 151 to evaluate higher-order 
diagrams involving double-excitations. These schemes can easily be generalized to 
handle triple-excitations and quadruple-excitations. 

Bubble Diagrams. — The bubble diagrams which are shown in Figures 13 and 14 
are required when the reference function for a closed-shell system is not defined by 
the matrix Hartree-Fock model, when a restricted Hartree-Fock reference 
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function is used for open-shell systems, or when molecular properties are to be 
evaluated. General programs have been written which allow any arbitrary one- 
electron operator matrix elements to be used in the evaluation of the bubble 
diagrams. 158 The expressions corresponding to each of the diagrams in Figure 13 
are as follows 

£ 2 (A) = £22 lv§S*L a » (121) 

y ab -^ijab 

£a(B) = 22 (122) 

i a -* la 

where A is the 'bubble' matrix. 

The expressions corresponding to the diagrams shown in Figure 14 are as 
follows 

E 3 (A) = i 2 2 (fi}abfcdabfc<ndl(.@ijab@ijcd) (123) 
ij abed 

E 3 (B) = - 2 2 (fi }a bfkaicfbc}k)l(&i]ab®ikbc) (124) 
ijk abc 

E 3 (C) = |2 E (fi}abfm]fabkl)l(®ilab@klab) (125) 
ijkl ab 

E 3 (D) = 22 (A ia A a1 ,A bi )l($ia®ib) (126) 

i ab 

E 3 (E) = - 2 X(AaiA iJ A i a)K® ia ®ia) (127) 
ij a 

E 3 (F) = 22 (A t aSibaiAy)K9ta®jb) (128) 
ij abc 

E 3 (G) = 22 (SiJabJ'at,c}A te )/(® tJ ab®t e ) (129) 
V oftc 

E 3 (H) =-22 (fijabj ~kbi]A k a)l(&i]ab®ka) (130) 
y'A aft 

£ 3 (I) = 22 (.SuabA ib A ta )K® t fal,@la) (131) 
(/ aft 

£ 3 (M) = -i 2 2 (fi]abA jk fi k ab)K2>ijab@ikab) (132) 
yfc aft 

^3(N) = i2 2 (J r ijabA b cfaci])K®ijab@i)ac) (133) 
y oftc 

The energy associated with diagram G is the complex conjugate of that 
associated with diagram J. Diagrams H and K, I and L are similarly related. 

Vector Processing Computers. — We conclude this section by considering briefly 
the impact which vector-processing computers are likely to have on many-body 
perturbative calculations in the near future. As noted in the previous subsection, 
use of the CRAY 1 has already enabled the triple-excitation component of the 
correlation energy to be evaluated for a molecule using a fairly large basis set. 188 

The algorithms devised to perform many-body perturbation calculations are 
well suited to vector processing computers such as the CRAY l. 141 > 158 The scalar 

158 S. Wilson and D. M. Silver, unpublished work. 

"» R. W. Hockney, Contemp.Phys. 1979, 20, 149; M. F. Guest and S. Wilson, 'Proceedings 
of the American Chemical Society Symposium on Supercomputers in Chemistry', Las 
Vegas, 1980. 
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operations on this machine are about twice as fast as those of the CDC 7600 or 
the IBM 370/195 whereas the measured time on the CRAY 1 for matrix multi- 
plication, using the vector order repertoire is twenty times faster than the best 
hand-coded routines on the CDC 7600 or IBM 370/195. 

The innermost loop in the program written to evaluate the triple-excitation 
energy component described in the previous subsection has the form 

DO 1 ID=1, N 

DD(ID) = D1(ID)*D2(ID) 

DE(ID) = D1(ID)*E2(ID) 

ED(ID) = E1(ID)*D2(ID) 

EE(ID) = E1(ID)*E2(ID) 

1 CONTINUE 

where Dl and D2 contain direct integrals and El and E2 exchange integrals. 
We performed a timing test for this loop and the results are shown in Table 2. 

Table 2 CPU times on the CRAY 1 and the IBM 360/195 computer required to 
execute the inner loop in the evaluation of the triple-excitation energy 
component 

Time b 

N a CRAY l c IBM 360/195* 

16 3.6 41.8 

32 5.2 82.9 

48 6.8 124.4 

64 8.4 165.6 

° DO-loop range defined in text. 6 CPU time in microseconds. c Using vector order repertoire. 
* FORTRAN H extended plus compiler with OPT = 2. 

The efficiency of algorithms based on diagrammatic perturbation theory 
together with the capabilities of the new generation of vector-processing com- 
puters will allow increasingly large basis sets to be employed over the next few 
years, thus reducing one of the major sources of error in present-day molecular 
calculations, the basis set truncation error. 159 

6 Some Applications 

General Remarks. — In this section we attempt to give a brief overview of the 
accuracy of electron correlation calculations which are possible at present using 
many-body perturbation theory. In the first part we describe some new calcula- 
tions on the lithium dimer and some previously reported work 92 on the nitrogen 
molecule. In the second part we review the application of many-body perturba- 
tion theory to potential energy curves. In part three, the importance of triple- 
excitations and quadruple-excitations is discussed in the light of recent calcula- 
tions. The calculation of molecular properties is discussed in the final part of this 
section. 
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Application to Li 2 , N 2 . — The results of calculations through third-order are 
shown in Table 3 and Table 4 for the lithium dimer and the nitrogen molecule, 
respectively. These calculations were performed by using the universal basis set 
containing functions with s, p, and d symmetry described in Section 3. For the 
lithium dimer the calculations were performed at an internuclear distance of 
5.0507 bohr, the experimental equilibrium value, whereas for the nitrogen 
molecule the calculations were performed for a nuclear separation of 2.0 bohr. 
The experimental equilibrium nuclear separation in the nitrogen molecule is 
2.068 bohr. The second-order and third-order energy components are given for 
both the model perturbation series and the perturbation series using shifted 
denominators. In addition to the simple perturbation results, [2/1] Pad6 approxi- 
mants and variational upper bounds are presented. 

The model perturbation series is observed to converge more rapidly than the 
shifted expansion, if we use the ratio of the third-order energy to the second- 
order energy as a measure of the rate of convergence. Note that for the total 
energies shown in Tables 3 and 4, the Pade approximants resulting from the two 
perturbation series are in closer agreement than the simple third-order energies. 
The difference between is [3/0] and £[2/1] is smaller for the model scheme than 
the shifted expansion. The is [2/1] values derived from the model scheme are thus 
of most interest although the optimized upper bound is also useful. The empirical 
correlation energy of the lithium dimer is estimated to be —0.126 hartree. For 
the model perturbation series the formation of [2/1] Pad6 approximants leads to 

Table 3 Components of the calculated correlation energy for the lithium dimer 
obtained using third-order perturbation theory and an s,p,d universal 
basis set a 

Model series Shifted series 



£ scf -14.871 23 -14.871 23 

£[2/0] - 14 . 969 20 - 14 . 993 76 

£[3/0] -14.98219 -14.983 74 

£[2/1] -14.98417 -14.984 49 

£var(A=l) -14.978 31 -14.968 99 

tfvarCAopt) - 14 . 979 20 - 14 . 972 76 

A 0P t 1.102 08 0.828 63 



a In hartrees. 

Table 4 Components of the calculated correlation energy for the nitrogen molecule 
obtained using third-order perturbation theory and an s,p,d universal 
basis set a 

Model series Shifted series 



£ sc{ -108.992 10 -108.992 10 

£[2/0] - 109 . 439 08 - 109 . 523 85 

£[3/0] -109.43676 -109.401 38 

£[2/0] - 109.436 77 - 109.424 31 

£var(A=l) -109.392 51 -109.335 28 

£var(A opt ) -109.396 54 -109.380 01 
A opt 0.904 8 0.729 5 



° In hartrees. 
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a calculated correlation energy of —0.1129 hartree which corresponds to 89.6% 
of the empirical estimate. Formation of an upper bound leads to a correlation 
energy of -0.1080 hartree which is 85.7% of the estimated total correlation 
energy. 

For the nitrogen molecule calculations were performed for nuclear separations 
of 1.75, 2.00, 2.25, and 2.50 bohr and these values interpolated to obtain a 
correlation energy at the equilibrium nuclear geometry. The empirical correlation 
energy for the nitrogen molecule is -0.538 hartree. The calculated correlation 
energy obtained by interpolation is -0.4501 hartree which represents 83.7% of 
the estimate. 

Both of the calculations described above recover a greater percentage of the 
correlation energy than previous studies. There are two deficiencies in the above 
calculations. The first, and probably the major effect, is that resulting from 
truncation of the basis set used in invoking the algebraic approximation. The 
concept of a universal even-tempered basis set should go a long way towards 
developing the large basis sets which are ultimately going to be required in order 
to achieve chemical accuracy (one millihartree). The basis sets used in the 
calculations described above do not include functions with <f> symmetry which 
will be required for an accurate description of angular correlation effects. 
The second deficiency results from the truncation of the perturbation series at 
third-order. Results of calculations for the water molecule will be discussed 
more fully in the fourth part of this section. We feel that it is always important 
to evaluate all terms, through whatever order the perturbation series is taken. 

Potential Energy Curves. — Potential energy curves have been calculated by many- 
body perturbation theory for a number of diatomic molecules, including FH, 160 
CO, 92 - 161 CS, 162 N 2 , 92 ' 163 BF, 92 F 2 , 163 Be 2 , 164 Mg 2 , 165 We take the results for the 
FH molecule as representative. 

The potential energy curves for FH are shown in Figure 15. Curves (d) and 
(e) in this Figure were determined by third-order many-body perturbation 
calculations. Full details of the calculations can be found elsewhere. 160 Curve 
(a) is the potential energy curve obtained by the self-consistent-field approxima- 
tion. Curves (b) and (c) were determined by Dunning 168 using the generalized 
valence bond method and by Meyer and Rosmus 167 using the coupled electron 
pair method, respectively. Curve (d) corresponds to the upper bound determined 
from third-order perturbation calculations using the model perturbation series 
and curve (e) to the [2/1] Pade approximant constructed from this series. Morse 
curves derived from experimental data are labelled (f) and (g). Curve (f) includes 
a correction for relativistic effects and represents an estimate of the nonrelativistic 
potential energy curve. 

«° S. Wilson, Mol. Phys., 1978, 35, 1. 

161 S. Wilson, Int. J. Quantum Chem., 1977, 12, 604. 

168 S. Wilson, /. Chem. Phys., 1977, 67, 4491. 

i«3 M. Urban and V. Kello, Mol Phys., 1979, 38, 1621. 

i«* R. J. Bartlett and G. D. Purvis, Int. J. Quantum Chem. Symp., 1978, 14, 561; M. R. A. 
Blomberg and P. E. M. Siegbahn, ibid., p. 583; M. A. Robb and S. Wilson, Mol. Phys., 
1980, 40, 1333. 

*" R. J. Bartlett and G. D. Purvis, /. Chem. Phys., 1978, 68, 2114. 

186 T. H. Dunning Jr., /. Chem. Phys., 1976, 65, 3854. 

167 W. Meyer and P. Rosmus, /. Chem. Phys., 1975, 63, 2356. 
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Figure 15 Potential energy curve for the hydrogen fluoride molecule 



The equilibrium nuclear distance determined from the self-consistent-field 
curve is 0.897 A which differs by 2.2% from the experimentally determined 
value of 0.917 A. Curve (e) yields an equilibrium distance of 0.908 A, which 
differs by 1.0% from the experimentally determined value. The fundamental 
frequencies of vibration derived from curves (a), (e), and (g) are 4471, 4261, and 
4137 cm -1 , respectively. 

Triple-excitations and Quadruple-excitations. — Only double excitations contri- 
bute to the correlation energy and other properties through third-order; triple-, 
quadruple and higher-order excitations are usually neglected. Triple- and 
quadruple-excitations first arise in the fourth-order of the perturbation series. 

To illustrate the importance of triple and quadruple excitations we present in 
Table 5 a full fourth-order calculation for the water molecule. 136 ' 138 - 166 This 
calculation corresponds to the equilibrium nuclear geometry and the basis set of 
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Table 5 Contribution of triple-excitations and quadruple-excitations to the cor- 
relation energy of the water molecule at its equilibrium nuclear geometry 





Triple 11 


Quadruple 




excitation 


excitation 


Diagram 


energy 


energy 


A 


-2.29 


-7.14 


B 


-1.33 1 


> -1.02 


C 


-1.35 J 


D 


-2.39 * 


> +3.51 


E 


-9.42 J 


F 


-7.11 1 


i- +7.86 


G 


-5.94 j 


H 


-7.77 




I(L) 


+ 2.36 




J(K) 


+ 9.73 




M(N) 


+ 1.53 




O(P) 


+ 1.15 




Total linked diagram energy 


-7.86 


+ 3.21 


Unlinked diagram energy 




-17.11 


E[2/l]* 


-0.2851 




E iS » 


-0.0020 






-0.0043 





• In millihartrees. 6 In hartrees. 

39 Slater functions first presented by Rosenberg and Shavitt. 168 The quadruple- 
excitation component can be divided into a part which corresponds to unlinked 
diagrams and a part which corresponds to linked diagrams. The unlinked diagram 
is only of interest since it represents the largest contribution to the size-con- 
sistency error in the widely used method of configuration interaction limited to 
single and double excitations with respect to a single-determinantal reference 
function. 

It can be seen that the linked triple excitation energy is in fact quite large and 
certainly chemically significant. There is clearly a danger in any technique which 
attempts to make a partial evaluation of higher order terms in the perturbation 
series. 131 ' 156 

Molecular Properties. — Kelly 5 showed that the many-body perturbation theory 
provides a powerful method for the calculation of atomic polarizabilities. 
Typical diagrams which arise in such calculations are shown in Figure 16. 
Itagaki and Saika 126 ' 127 have extended Kelly's work to molecules. The present 
author 120 has demonstrated how the use of a scaling parameter can lead to 
improved results in such calculations. Itagaki and Saika have also very recently 
reported applications to nuclear spin-spin coupling constants. 169 

Many-body perturbation theory can also form the basis of a useful tech- 
nique for the calculation of photoionization effects in atom and molecules. 
Kelly has recently reviewed this area of research. 170 



168 B. J. Rosenberg and I. Shavitt, /. Chem. Phys., 1975, 63, 2163. 
189 T. Itagaki and A. Saika, /. Chem. Phys., 1979, 70, 2378. 
170 H. P. Kelly, Comput. Phys. Comm., 1979, 17, 99. 
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Figure 16 Diagrams which arise in the calculation of polar liabilities. The heavy dot 
represents a dipole moment matrix element 



General Remarks. — In this final section, some other aspects of the many-body 
perturbation theory of molecules are briefly reviewed. The final comments 
address the general properties which we believe methods for treating electron 
correlation in molecules must possess. 

Some Other Aspects. — In 1972, Coulson (quoted in ref. 6) suggested that dia- 
grammatic techniques would become increasingly important in theoretical chem- 
istry. The impact of diagram methods is perhaps only just being fully realized. 
In this section very recent work and possible directions for future research are 
discussed. In particular an extension of the work described in the second part of 
section 4 on invariants is described, the possibility of employing group-theoretical 
developments in many-body perturbation theory is indicated, recent work on 
electron-molecule scattering is mentioned, and the calculation of energy deriva- 
tives is briefly discussed. 

Even-order Invariants for Rayleigh-Schrddinger Perturbation Theory. Full fourth- 
order correlation energy calculations are now possible and it is, therefore, 
important to have fourth-order forms which share the invariance properties 
of the [2/1] Pade approximants in third-order. The expression 



7 Concluding Remarks 




(134) 




171 



(135) 



171 S. Wilson, Int. J. Quantum Chem., in press. 
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Graphical Methods of Spin Algebras. In most problems of chemical interest, 
the molecular hamiltonian can be taken to be spin-free and can be written 172 ' 173 
in terms of the generators of the unitary group, U(n), where n is the number of 
basis functions 

J? = 2 <a\z\by + \ 2 <ab\v\cd> (#ac*bd-h>c#ad) (136) 

ab abed 

in which 

#ab = T,aZ a a ba (137) 

a 

and a + and a are the usual creation and annihilation operators. A spin-adapted 
perturbation expansion can be obtained by defining 

= S<flt|/|fl>*«. (138) 

a 

and 

l#\ = i L <ab\v\cd>(<?ac<?bd-6bc#ad) + S {<a\z\by^ab-Sab<a\f\b}S' ab } 

abed ab (139) 

The matrix elements of the generators, & &, can be evaluated using the graphical 
methods of spin algebra. 174-176 The particle-hole formalism introduced by Flores 
and Moshinsky 177 and recently discussed by Paldus and Boyle 178 is most suitable 
for such developments. 

Electron-Molecule Scattering. Kaldor and co-workers 179 have recently shown 
that diagrammatic perturbation theory can be used to expand the optical potential 
in electron-molecule scattering calculations. The second-order diagrams, in- 
cluding those resulting from electron exchange are shown in Figure 17. 

Pade Approximants to Energy Derivatives. The evaluation of derivatives of the 
energy, with respect to the nuclear geometry, forms the basis of an efficient 
approach to the calculation of potential energy curves and surfaces. It is easily 
shown that 

$* - ^ J 2 (iZl) <M- 1)""* £ &<P) (140) 

where 

P = (141) 

172 J. Paldus, Theor. Chem., Advances and Perspectiv es, 1976, 2, 131. 

173 I. Shavitt, Int. J. Quantum Chem., Symp., 1978, 12, 5; and in 'Proceedings of Daresbury 
Study Weekend', November 1979, ed. M. F. Guest and S. Wilson, Science Research 
Council, 1980. 

174 A. P. Yutsis, I. Levinson, and V. Vanagas, 'Mathematical Approaches of the Theory of 
Angular Momentum' Israel Program for Scientific Translation, 1962. 

175 E. El Baz and B. Castel, 'Graphical Methods of Spin Algebras', Dekker, New York, 1972. 

178 S. Wilson, in 'Proceedings of Daresbury Study Weekend', November 1979, ed. M. F. 
Guest and S. Wilson, Science Research Council, 1980. 

177 J. Flores and M. Moshinsky, Nucl. Phys., 1967, A93, 81. 

171 J. Paldus and M. Boyle, Physica Scripta, 1980, 21, 295; M. Boyle, Dissertation, University 
of Waterloo, 1979; J. Paldus, in 'Proceedings of Daresbury Study Weekend', November 
1979, ed. M. F. Guest and S. Wilson, Science Research Council, 1980. 

179 A. Klonover and U. Kaldor, J. Phys. B: Atom. Mol Phys., 1978, 11, 1623; ibid., 1979, 12, 
323. 
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Figure 17 Second-order diagrams for electron molecule scattering 

Ek is the kih order energy corresponding to the zero-order hamiltonian jfi and 
Etc is the kth order energy corresponding to /"^o- It therefore follows that the 
[N+l/N] Pade approximants to the energy derivatives share the invariance 
properties of the [N+ l/N] Pade approximants to the energy itself. 

Final Comments. — Many-body perturbation theory has several theoretical pro- 
perties that make its use in molecular studies attractive. If a system is composed 
of non-interacting subsystems, then this fact should be evident at any level of 
the theory. It is possible to formulate any quantum mechanical problem in such 
a way that, if the whole system is separated into subsytems, all quantities appear- 
ing in this formalism are additive. There is a close connection between this 
question of additivity of properties and linked diagram expansions and Lie 
algebraic formulations of quantum mechanics. Many-body perturbation theory 
provides a formulation of the correlation problem which satisfies these require- 
ments. The energy of any system is always written as a sum of the energies of its 
component parts. The energy is written in terms of single-particle state functions 
rather than iV-electron functions. These theoretical properties lead to tractable 
computational schemes for treating correlation effects with little more effort than 
a conventional self-consistent field calculation. 180 

Coulson 181 thought that it would take 15 years for the impact of diagrammatic 
techniques to be fully realized in theoretical chemistry. It is therefore not surpris- 
ing that the first half of this period has been devoted mainly to the development 
of new methods and algorithms. Although this development will undoubtedly 
continue, it is clear that, armed with these new techniques, theoretical chemists 
will be able to attack problems with an accuracy which was not previously 
attainable. For example, it will be possible to calculate rotation constants for 
small molecules more precisely. Although the accuracy of calculated rotation 

180 m. F. Guest and S. Wilson, 'Proceedings of the American Chemical Society Symposium 
on Supercomputers in Chemistry', Las Vegas, 1980. 

181 C. A. Coulson, 1973, Progress Report, Department of Theoretical Chemistry, Oxford, 
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constants will still not match that of microwave experiments, calculations are of 
great value when experiments are difficult or impossible. In the past few years, 
theoretical studies have been of great value in the detection of certain interstellar 
radicals and ions 182 that are very difficult to study in terrestrial experiments. 



™* S. Wilson, Chem. Rev., 1980, 80, 263. 



